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MULTISPECIES GOMPERTZ-TYPE MODEL WITH TIME DEPENDENT 
PARAMETERS AND ITS STOCHASTIC GENERALIZATION 


S. S. DE 


Department of Applied Mathematics Calcutta University, University College of 
Science, 92, Acharya Prafulla Chandra Road, Calcutta 700009 


(Received 12 July 1985; after revision 3 April 1986) 


A deterministic model for several-species-population which follows Gompertz’ 
law, has been proposed and extended for the case of time-dependent periodic 
and non-periodic growth and interaction rates. The method of solutions has 
been given together with the explicit solutions for the two-species population. 
The stochastic generalization of this model has been proposed. The Fokker- 
Planck equation corresponding to the stochastic differential equations arising 
out of this model, has been approximately evaluated for the transition pro- 
babilities for both the stationary and non stationary cases. 


1. INTRODUCTION 


The deterministic models of population growth and of the interaction among the 
biological species are based on the phenomenological equations such as Verhulst 
Logistic equation, Lotka-Volterra equations and Gompertz’ Law*’®. Again, in the 
stochastic model the environmental fluctuations on the macroscopic behaviour of the 
linear and non-linear system have been considered. Several mathematical tools are 
employed there; for example, the Markov process techniques and the stochastic dif- 
ferential equations, together with their conversions to the Fokker-Planck equations. In 
some cases, the interaction among the species has been taken as the interaction of 
Markov processes, whereas, in some other cases, the ‘stochasticity’ has been introduced 
by inserting the random parts in the parameters of the deterministic models 
and thus the stochastic differential equations are generated. In fact, in an 
earlier paper®, the stochastic version of non-linear model of several interacting species 
in a randomly varying environment together with predator-prey type random interac- 
tion, has been investigated. There the Fokker-Planck equation arising out of the 
stochastic differential equations has been integrated in the form of the path-integral 
which in turn has been approximately evaluated for the transition probability for the 


stationary case. 


In a previous paper’, we have generalized a deterministic linear model for 
n-species biological population, originally proposed by Coutlee ef al.” and later by 
Gomatam® for two species interaction. There we have proposed the interactions 
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among the species, which follow Gompertz’ law and discussed the situations where 
such models may be relevant. Also, the stochastic version of the model through the 
introduction of white noise type random components in the parameters of model has 
been investigated. The transition probability has been evaluated approximately under 


the steady state. 


In the present paper we consider the generalized model, proposed in the above 
mentioned paper, for the situation when the growth and the interaction rates become 
time-dependent. Also we form the stochastic models from these deterministic equations 
and investigate the resulting stochastic differential equations. The corresponding Fokker- 
Planck equation has been solved approximately for the transition probability for the 
stationary case. In section 2, we introduce and investigate the model when some or 
all of the growth and the interaction rates are time-dependent. In the subsequent 
section 3, the stochastic model arising out of this deterministic model through the in- 
troduction of random comp onents in the parameters, has been studied. 


2. THE DETERMINISTIC MODEL WITH TIME-DEPENDENT RATES 


The governing equations of the Gompertz-type-interaction model, which was 
proposed by the present author* in the previous paper to describe the interaction among 
the several biological species, are 


ou = 2 bijy; + Zi for t => 0 } 
Zo ME mi) 
where , 
yw =nX (Fes 19 20 en) | 
and 
vs = y (0) at t = 0, for all i ...(la) 


Ay 8 (bee le 2 1 hes n) and b1,’s (i, j = 1, 2, ..., n) 


are the growth and the interaction rates re 1 i 
| spectively. X; is the populati i 
the i th species at time tf. e en 
We now consider this deterministic 


: ’ population model with time- 
and interaction rates. In this case, the Reco PaRe cmb erow ty 


governing equations are 


Q 


d= Py Oy +20) 0) 


(i = 1, 2, ..., 2) 


with the initial conditions 


yi =y(0) att =0. 
...(2a) 
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Now let 5;; (t) and Z; (t) be periodic functions of time with period +. Such periodic 
nature of the growth and interaction rates may arise out of the seasonal changes in 
the environment or from periodic nature of mating of animals and subsequently the 
breeding of those species. Also, such fluctuations in the interaction terms occur due 
to periodic variations of common resources such as food, etc. Then setting w;(t; +) 
= yi (t + +) — yr (t); we have the following equation for w; (t; t): 


A oe . by; (t) w; (t; 1). -os(3) 


Let us specify the periodic functions b;, (t) and Z; (t) as, 


bi; (t) = (1 + € cos pt) bi; 
and 4) 
Zi (t) = (1 + € cos pt) Z: 
with p = 2n/t. b,,s and Z, ’s are independent of time and € is small. 
Then we have from (2), 
(1 + € cos pt-’) oo = % bi; w; Sa BY 
where 
B= site ; 
with Zi = = by %; -.-(6) 
j 


Here «,’s can be thought as the equilibrium values of », ’s for a model with time-in- 
dependent growth and interaction rates. Using (4) we have from (3), 


| dw, (t; 7) 


cial wet is belie). 0 
1 + € cos pt dt ; err ( 


The two set of equations (5) and (7) are of the same type. The solutions of them [for 
(5)] can be obtained by substituting 


Pee is oa ...(8) 
Pp 
Then eqns. (5) becomes, 
dw ca 
ap = 2 ns 9) 


The solution of these equations are of the form : 


B= SAV, ...(10) 
t 


1084 S. S. DE 


where V! are the eigenfunctions of the matrix {b,,} with the eigenvalues A; and A; are 
‘ g j 


constants with respect to T. In the case of multiple roots of the characteristic or secular 
equation for the matrix {b;,}, (that is, the case of some or all of the A} s are equal, 
the above solutions are somewhat involved. In fact, in this case, the so-called similar 
terms appear and the constant eigenvectors vi , in the formula, are to replaced by a 


polymomial in 7, which are of the type 
Vi tet... 
Also, it should be noted that if the matrix {b; js are symmetric, then the eigenvalues 


A; are all real, but in the general case A; may be complex. Now the solution becomes 
(for the case of A, ’s are all the different) 


wy = > Al Vi exp (2 (: | ae )) mile) 
F 


or 


Mz=m—ayt+ Ss A; V; exp ey (: +e ar Yh. ...(1 1a) 
t 


The constants A, ’s can be determined from the initial conditions (2a). 


In the case when the interaction rates b 4; S are constants whereas the growth 


tates Z; (t)’s are perodic functions of time with period t, we have from (2) 

dy, 

LI elon 2 bi; Yj + Zy (t). : avel ta) 
Integrating, from ¢ to t + +, we have, 


yr +r 
P t is = , , 
y(t + 2 y(t) = bi; Sa (ty’) dt’ + J Ze (t’) dt’. cfCTS) 


Now, using the solutions of (3) when b yy 'S are independent of time 


Mit} t) = w(t+)—y(t) = 5B PV! e 
1 t 
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we have, from (13), 


1 . t45 tyr 
EBV, eM=Eb, fy (r)dt + f Z(e)de 
t 
From this it follows that, 


far 
I 1 
=| Ye (t') dt’ = a 5 By (b> a V eA! 
Le 


t 


t+r 


= 3 (O")is +| Z, (t') dt’ (k = 1, 2, ..., nm) ...(14) 


where (b~"),: are the elements of the inverse matrix of the matrix {b,,}. 


These equations give the time averages over the period t of the population sizes. When 
for all the eigenvalues A. Re A; < 0, we have, as t — oo 


<—Yo>r = — - (bai <Zi>r 


(hie Ty, 7) see 29) 
where < >, denotes the average over the period t. 


Now we proceed to solve (12) for the case when Z; (t) are arbitrary functions of 
time. Following Gantmacher*, we use such a linear nonsingular transformation of 
variables 


y=U0 (det U # 0) ..(16) 


that, U~! BU reduces to the “triangular” form, 


f Ay Ci. Cis : Cin) 
| 0 Ao Cys : Con | 
Bw 0 As: Ce, | 
U2 BU = | j | = C(say) ati) 
Wee ee Mea tS, 
Wicca lit Aaa ie Pee 


where A1, As, ..-, A, are the eigenvalues of the matrix {b;;} = B, while the moduli of the 
nondiagonal elements Cy. (i < k) may be made arbitrarily small by suitable choice of 
transformation (16). 
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By using this transformation, equation (12) reduces to 


Shan Z (1) 
where f ...(18) 
Z(t) = U Z(t). / 


The last equation of this matrix equation is 





dz, ~ 
dt = An %m + Z, (t). ...(19) 
Its solution is 
; a } 
dy (t) = e4,' § e-An! Z, (t) dt + eAnt 7,(0) 
0 
. ...(20) 
where | 
v(0) = U- ¥(0). J 


Substituting 2, (1) from (20) into the last but one equation of (18), we get after solving 
for 2,-1 (t), 


t ~ 
Dy-1 (t) = eXn-1' f e7*n-a" Z,_, (t) + eAn-2! ¥,1(0) 
0 


i t t ~ 
+ Canton Obnat J Anat Ant {I e-Ané Zy (1) dt + U,(O)}dt. 
eQn 


In this way we can solve, exactly, all the equations in (18) for 2,_, (f), Deca (t), oo: 


My (t). 


We now like to find the explicit solutions for a two-species population with time- 
dependent periodic and non-periodic growth rates. We take some specific values for 
the interaction coefficients, b}; = —1, b., = —1, b,, = 0.09 and b.. = 0. The value 
of b,. = 0 corresponds to the zero value of the self-limiting term of the Som speci 
which helps to the growth of the first; whereas the first species is hostile to the = eF 
species. The U-matrix and its inverse, for this case, can be found to be Bare 

> 


0.1 — 
u =| 5) and u = ( 0 et) 
1/0 yen P04 -» (22) 
(€ Small) 


Then matrix C becomes 


bw i — 0.1 € 
te ad be ( mathe 


0 —0.9 
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which is a triangular matrix with small off-diagonal element. Then eqns. (18) become 


; $ Z, (t) 
(2 )- (06'S) + (Sx@sar20 ) 20 


Ve 
a, 
> ae eet + 0.1 ye : 
% : 


We find the solutions of these equations for the first case of periodic growth rates : 
Zi (t) = pil + €, COS pt) (i = Lap. we) 


The solutions are 


where 





Vi (t) = 0.1 e's Yo (0) + gL (e-orl == e709") 


+ fx (1 — e714) 


; (p sin pt + 0.1 cos pt) we é 
1 ps 0.1 + 10p* sas ee 


(1 ps oun" Ir) Pi Capen — e- ol) 
x{ 0.9 hae \ 
0.9 Cty = e-o-1t) 
8 (0.81 + p*) 
p sin pt + (0.09 — p?) cos pt } 
+—~0.81 + p*) (0.1 + 105?) 


Tt pi €é, + 0.1 fe «) | 





a * eae! 
+ (fi — 0.1f21) —p.9-— + (pier — O.1p, <,) 


Ake sin pt + 0.9 cos pt) — 0.9e-°" 


...(26a 
0.81 + p? (26a) 


y(t) = e973 (0) + BOLO MO (044 — g-0my 


= sin pt + 0.1 cos pt) 
+ 10f, (1 — e-°™) + poe ws 0.01 + p? 


nee O°-9f __ e-” yi 
“8 (0.81 + p’) 


(p sin pt + (0.09 — p*) cos pt } 
(0.81 + p*) (0.01 + P*) 


(equation continued on p. 1088) 


+ (— pat 0.1p, a) {= 


a 
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1 J 
roe jr +017) 4 O99 (ect) a ne 


— e-01t) } ...(26b) 


—0-% 


For the second case of non-periodic growth rate (a damping oscillatory) of the type 
Z=piteae* (i= 1, 2) ued) 
with a having posilive real part, the solutions become | 


y1 (t) = 0.1 yo (0) e- ole 4 ee | (0) . O.1y. (0) (e704 ak) e~0°9) 





(= hr + O.1fs) a ea ont Cre eee 
a rae 0.1 0.8 


(— «, + 0.1 €) { enn e-0m } 


2oeae si a VAR oe 


+ peo (hee. 2!) 
—at 
+ 2 Fog + Or (0) — 0.1y (0)) 2-0 
«, — Ole, 


Tat OO aa tac) eae ...(28a) 


Ya (t) = 20) ee + —HO FOI (eo _ gon) 





4 (=p; + 0.1p,) LO 9-0-8 _ gros 
0.9 








0.1 ios i ee 

(heater meer on 

410%, (Le ett ee, rae ae ...(28b) 
For large t, the behaviour of these solutions are, for the first case, 

mim pe {ite PSDB os pe | 
+p eines {091 Asngs = repr 
Y(t) = Sy Oe — 3) + e's eer eee 
pe +01 ja) {PABLO peogeh 


...(29) 
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yn) = pe 1 
1007. Se ...(30) 
yve(t) = heQe (po — pr). ! ( 
J 
Thus we see that, in the first case, the asymptotic behaviour of the soluations is peri- 
odic about the equilibrium solutions of the case with the time-independent parameters; 


whereas for the second case the solutions reach the equilibrium values asymptotically, 
that is, to the asymptotic stable values. 


3. DISCUSSION ON STOCHASTIC GENERALIZATION 


In a previous paper’, we have already made the stochastic generalization of the 
model with constant growth and interaction rates. Now we intend to generalize the 
present model by assuming that the interaction rates have random parts. Thatis we 
can write the governing equations as, 


dy r 
a = * (G1; + cry 7 (O)) Yi + Zi (0) - (31) 


where y, (t) are 5-correlated Gaussian white noises, 


with 
<y,;()> =0 
and ...(3la) 
<yi(t +) yj (t)> = 51; 3 (*) 
and 6,,’s are the averages of b,,’s. We may write (31) as follows : 
: a 
yi = Ki(y) + = 81; (Y) x; (0) | 
with | 
Ki (vy) = bij y, + Z(t) : ...(32) 
j 
| 
and | 
Bij (Y) = Cry; J 
The corresponding Fokker-Planck equation can be written as** 
o 1 D3 oe 
=- —_ — ) — Sa Saeed B v P eee 33) 
Be. jy, de PI + 5 By, Byy Bo» (y) PI ( 
& Bey 
with 
Ap) = BOpy +4 Cuy Cy) Vy + Zp (0) 
and .-.(34) 


—-—~-— — SY 


Buy (Y) == Cu; Cy Y5. 
J 
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The path integral solution of (33) can be found and the approximate expression for the 
the transition probability P (yi, T | ¥:,) can be written following the same method of 
the previous papers**. They are, respectively, 


00 ti 7 
P(y1, Ti | 0) = {det BOAT)? J D (y)exp{— 3S L dt 
ae | 
where f . t 
L=[y— AW BOYD — AQN+22 Cy | ...35) 
| 
AQ)=E Cy — 5 Cy + BO 
and 


P (yi, T | yi, 9) = | det B (0° (T) ee 


exp \- a L(y, 2), ) at } j DG) exp{— 44 0 (&) dt} 


0 


...(36) 


where y~, ’s are the solutions of the equations 


yr — * qij Vj 4+ Zi (t), (i = L ae eeey n) 
with f pa 
qij = bij _ 3 Ci; Cj; G7 == ie yy oeey n). J 


The solutions of these equations can be obtained in the same way as for the equa- 
tion (12) of the deterministic case and they are of the form (26) or (28) for the two- 
species populations, as for the example. From (36), we can then write down the ex- 
pression for the transition probability using these solutions y; if ¥ 


For a particular case, when Cy; = ¢,; + C 5,,; with ¢,,;’s being very small com- 
pared to C; we find, 


qij = bi; — B{C* Sy + Ci; €j, + 41} ) 
and 


Buy (y) = C? Buy y?. .- (38) 


-—+— 


In this case the transition probability becomes 


P (yi, c& | y1, 0) = D exp \- o In| det B (> ()) }} ...(39) 
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where D includes the path integral part for &’s together with a constant term from 
c . : . 
L(¢ VeL, y ) , less important for the present discussions, 


We have, using (38), 


P (i, 1990) = D ear Il 





ex =) ry ...(40) 


When Z (ft) is a periodic function with small amplitude, yé (©e)’s do not attain 


fixed values. On the contrary they oscillate with small amplitude around some fixed 
values. If Z(t) is some slowly varying damping oscillatory function of time, yn ’S 


can attain asymptotically stationary values. In this case the transition probability 
has some definite value. 


In the present paper we have studied the deterministic models with time-depen- 
dent parameters. The methods of solutions have been discussed and the explicit solu- 
tions are also given for the two-species biological population. These models may also 
be relevant to other actual situations, more precisely to some ‘Controlled’ situations; 
as they are found to be relevant, for the cases of the time-independent parameters. 
discussed in the previous paper’. This may be investigated in future. 
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NOTE ON A THEOREM OF ALTMAN 


B. Watson, B. A. MEADE AND C. W. Norris 


Department of Mathematics and Statistics, Memorial University, St. John’s 
NF, Canada A\C 5S7 


(Received 6 January 1986) 


The purpose of the present study isto clarify and Generalize a theorem of 
Altman’. 


Let Q be a real-valued function satisfying 
(a) 0 < Q(s) < sfors < Oand Q (0) = 0, 
(b) g(s): = s/(s — Q(s)) is nonincreasing on (0, co), 


4 


(c) f g (s) ds < co for each s, > 0, and 
0 


(d) Q(s) is nondecreasing. 


We suppose throughout that Y is a complete metric space. A function F: X¥ ~X 
is called a generalized contraction if d (Fx, Fy) < Q (d (x, y)) for all x, yEX. Altman’ 


has shown that functions Q, having all four properties above, arise naturally in a 
Banach space setting. 


The purpose of this note is to clarify and generalize the following theorem of 
Altman’: 


Theorem 1—A generalized contraction has a unique fixed point z and, for every 
XE X, lim (F)" (Xo) = Zz. 


n->oo 


Actually, Altman’s version has the following condition on Q replacing (a). 
(7)0 < QO(s)< sforO<s<s. 


This form does not give uniqueness since s,; depends on the starting point Xp. 


However, (a) is sufficient for uniqueness and the rest of the proof goes through as in 
Altman’. 


We generalize Theorem | for a pair of mappings, one of which is a generalized 
contraction and the other expansive. 


Theorem 2—Let F, G: X > X be such that 
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(1) Fis a generalized contraction, 
(2) d(Gx, Gy) > d(x, 5) for all x, y € X, and 
(3) F(X) C G(X). 
Then F (x) = G (x) has a unique solution z and, for every X» € X, 
lim (G" F)" (x,) = z. 
nc 
PRooF : First, G~ exists since (2) implies G is one-to-one. Also, (3) gives 
G'F:X—>X. Next it follows from (2) that 
d(GG" x, GG! ») > d(G" x, G-'y) 
or 
d (x, y) > d(G"" x, Gy), 
Replacing x by Fx and y by Fy in the latter inequality gives 
d(G" Fx, G Fy) < d(Fx, Fy) < Q (d(x, y)). 


Hence G~ F is a generalized contraction with respect to Q. The result now follows by 
Theorem 1. . 


We note that taking Gto be the identity mapin Theorem 2 gives Altman’s 
result. 


In order to illustrate Theorem 2, we give the following example. Let ¥ = [0, 00) 
and let d be the absolute value metric. Define F,G : [0, oo) > [0, co) by F (x) 








= 1 pees Kc a SaaS ° <8 

volt get er Bike ok Let () $s) = 7416 Clearly Q satisfies conditions (a), 
yes oe ok ise cel tal 

(b), (c) and (d). Then d (Fx, Fy) = xy + a(x + y) + 16 S Se Ree 


= QO (d(x, y)) and d (Gx, Gy) =2|x—yl >|lx-—vyl. Also rx=(0, = | 
C [0, co) = GX. 


] . 
Hence, all conditions of Theorem 2 are fulfilled and 2x = xq 4 has a unique 


solution in [0, oo). 
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SOME RESULTS IN THE FIXED POINT THEORY, IIl 


MILAN R. TASKOVIC 


Prirodno-matematicki fakultet, Studentski trg 16, 
11000 Beograd, p. box 550 Yugoslavia 


(Received 20 December 1985) 


In this paper we give applications of our Localization Monotone Principle 
in the fixed point theory. 


Let X be a topological space, T: X > X and let B: X > R%: =[0, co) bea T- 


orbitally lower semicontinuous function on X. A function B mapping X into the 
reals is T-orbitally lower semicontinuous at p if {x,} is a sequence in O(x): 
= {x, Tx, T? x,...},x € X and x, > pimplies that B(p) < lim inf B(x,). A space 
X is said to be 7-orbitally complete if every Cauchy sequence which is contained in 


O (x) for some x € X converges in X (c. f. Taskovic"’). 


In connection with this, in a recent paper’* we introduced the concept of TCS- 
convergence in a space X; i. e., a topological space X satisfies the condition of TCS- 
convergence if there exists a point x € X such that B(T" x) > 0 (n > co) implies 
{T” x} has a convergent subsequence. We stated a new fixed point principle : 


Localization Monotone Principle'’—Let T be a mapping of a topological space 
X into itself, where X satisfies the condition of TCS-convergence. Suppose that there 


exists y € XY and a mapping 9: Ry — Rj such that (1) < rand lim sup 9 (z) < 1 
: zZ>t+0 
for every ¢ > 0 and 
B (Tx) < ¢ (B(x)), for every x € O(y) (1) 


where B: X > R is T-orbitally lower semicontinuous a i i 
Then 7 has a unique fixed point — € Y. ay Ae 
We get the following from the preceding principle : 


Corollary 1\—Let T be a mapping of a to 


ological into i 
that there exist « € [0, 1) and Pollogica! space X into itself. Suppose 


B (Tx) S « B(x), for allx € xe 
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where B: X¥ — R& is T- -orbitally lower semicontinuous and B (x) = Oimplies Tx =x. 


If for some x € X the sequence {7" x} has a convergent subsequence, then T has a 
fixed point & € YX. 


Corollary 2—Let T be a mapping of metric space (X, P) into itself and let Y be 
T-orbitally complete. Suppose that there exist « E€ [0, 1) and 
A (Tx, Ty) € « A (x, y), for all xXxVEX 


where A: X¥ x X¥ > R{, x >A (x, Tx) is lower semicontinuous and P [x, y] 
< A (x, y) for all x, y © X. Then T has a unique fixed point § € X and 7" x > — for 
eachx € X, 

We notice, this result is generalization of the results of Chakrabarty’. 

In this paper, we apply our Localization Monotone Principle in the fixed point 
theory. 

In this part we assume that (YX, P) is a complete metric space and that ® isa 
mapping of R° into itself satisfying the conditions: ® is nondecreasing and conti- 
nuous on the right, © (¢) < t for allt > 0. Further we assume that ¥ is a continuous 


mapping of ¥ x X into R$ satisfying the conditions Y (x, x) = 0 and P (x,y) 
< VY (x, y), for any x,y in X. 


Corollary 3 (Chakrabarty?)—Let M be a subset of Y and T be a mapping of M 
into itself. Suppose that 


¥ (Tx, Ty) < © (¥ (x, y)), for allx, y E M cet ay 


where ¥ (x, T* x) < A(x), € N for every x € M and A (x)a positive number. 
Then there isa point & in ¥ such that for anyx € X¥,T"x >~&(n > 00). If Mis 
a closed, then & is the unique fixed point of T. 


PrRoor: Let B(x) = ¥ (x, Tx) and p(t) = O(f). It is easy to see that B and 
g satisfy all the required hypotheses in localization monotone principle. Also, from (2) 
X satisfies the condition of TCS-convergence, because Y is a complete metric space. 
This completes the proof. 


Popa’ obtained a fixed point theorem which is a common generalization of 
results of Banach, Ray’® and Jaggi’. In this paper we extend this statement. Namely, 
we give applications of localization monotone principle, and as special case results 
of Popa’, also and some others results. 


Corollary 4 (Popa*)—Let T be a continuous mapping of a Hausdorff spaces X into 
itself and let f be a continuous mapping of Y¥ x X into the nonnegative reals such that 
for some a,b € Ri, a + 6 < | and forall x,y € Xandx + y: 


( (Tx, Ty) & af (x, Tx) f(y, Ty) ( f (x, y))-? + OF (x, Y) =e) 
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and 
f(x OLY) <f? (, y) and f (%, y) # 9. (4) 


If for some x € X the sequence of iterates {T" x} has a convergent subsequence, then 
T has a unique fixed point. 


Proor: Let x be an arbitrary point in ¥. Then, for y = Tx ~ x, from (3) 
we have 


f (Tx, T?x) < (a + b) f (x, Tx), for all x € p « 


Let B (x) = f(x, Tx) which is lower semicontinuous on X, and let 9 (t) =(a+5) t, 
for t€ R,; then B and 9 satisfy all the required hypotheses in our localization mono- 
tone principle. Since X satisfies the condition of TCS-convergence (the sequence of 
iterates {7" x} has a convergent subsequence), applying our localization monotone 
principle we obtain TE = § for some & € X. Uniqueness follows immediately from 
conditions (3) and (4). 


Corollary 5 (Jaggi’, Bose and Mukherjee*)—Let T be a mapping of a metric space 
X into itself satisfying the following condition 


P [Tx, Ty] & a P [x, Tx] P [y, Ty] (P (x, yD? + 5 P Ix, yl, x AY 


where a and b are nonnegative real numbers such that a + b < 1, If there exists an 
orbit O (x,) which contains a convergent subsequence of which T is orbitally conti- 
nuous, then JT has a unique fixed point. 


PROOF : is analogous to the proof of the preceding statement of Popa’. 


We note, this statement of Bose and Mukherjee’ is special case of Corollay 3 of 
Popa’. 


Let N denote the set of all nonnegative integers. A pair (X, —) of a set X and 


IN 
a subset ~ of the set X x Xis called an L-space if the following two conditions 
are satisfied : 


(i) Ifx, = x € X for alln E N, then ({x,}, x) « >. 
(ii) If (s, x) « +, then (t, x) « — for every subsequence tf of s. 


In what follows, we shall write s + x or x, -> x instead of (s, x) ¢ —, and read s 
converges to x, where s = {x,},n © N. If s = {x,},n © N,isa sequence in a set X, 
and if f is a mapping On X, then f(s) denotes the sequence { f(x,)}. Let (X, +) be - 
L-space. It is said to be separated if each sequence in XY converges to most one 
point of X. Let d be a nonnegative extended real valued function on X x X. The 


L-space (X, —) is said t - i i i 
: pace ( ) issaid to be d-complete if each sequence nen in X with 


i d (Xn-1, X,) < e© converges to at least one point of X. 


n= 
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The following statement is corollary of our localization monotone principle: 


Corollary 6 (Kasahara*)—Let (X, >) be a separated L-space which is d-complete 
for a nonnegative extended real valued function d on X x YX. , and f be a continuous 
mapping of X into itself satisfying the following conditions for some a, B with 
O<a<landO<B< oo: 


d (f (x), f2 (x)) < ad (x, f (x)) < « B, for allx € X. ..(5) 


Then f has a fixed point, and the sequence {7x} converges to a fixed point of f. 


Proor : As can readily be seen by induction, d (f" (a), f"*! (a) < «" d (a, f (a), 
for every n € WN, where N denote the set of all nonnegative integers. Hence the 
d-completeness of the space implies that the sequence { f" x} converges to some § € X. 
This implies that X satisfies the condition of TCS-convergence. Let B (x) = d (x, f (x)) 


and 9 (t) = at. Applying localization monotone principle we get § = f& for some 
EFEX. 


Also, an immediate corollary of the our preceding localization monotone 
principle is : 


Corollary 7—Let T be a mapping of metric space (X, P) into itself. Suppose that 
there exists « € [0, 1) such that 


P [Tx, T*x] < ap [x, Tx], for allx € X 


where x — Pp [x, Tx] is T-orbitally lower semicontinuous and {T" x} has a convergent 
subsequence. Then T has a fixed point in X. 


Corollary 8 (Dhage*)—Let T: X > X bean orbitally continuous self-map of a 
metric space X and let X be T-orbitally complete. If T satisfies the condition 
min {p [Tx, Ty], P [x, Tx], P Ly, Ty]} 

+ amin {p [x, Ty], P [v, TX], < p P [x, y] + 9P [x, Tx] ...(6) 
for all x, y © X anda, pandq are real numbers such that 0 < p + q < |, then for 
each x € X, then sequence {7” x} converges to a fixed point of T. 

Proor: Let x € X be an arbitrary point in ¥. Then, for »y = Tx, from (6) 
we have 


min {P [Tx, T?x], P [x, Tx]} < (p + q) P [x, Tx]. 


Hence, for B (x) = P [x, Tx], 9 (t) = (p + 9) 4, and since X satisfies the condition of 
TCS-convergence (X is T-orbitally complete metric space and P (Px, Pon's] < (p+q)" 
(1 — p — q)" Pp [x, Tx]), applying our localization monotone principle we obtain 


TE = & for some & € X. 
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Some Remarks 


(1) We notice, that and some other results (Theorems 2, 3 and 4) of Dhage’, are 
special case of our localization monotone principle. Proofs are analogous of the 
preceding proof of Corollary 4. 


(2) Also, in a paper'® we have proved the following statement, which generalizes 
the preceding statement of Dhage*. A comparative study of these generalization has 
been made more recently by Yeh?® and Khan’. 


Corollary 9 (Taskovic'*)—Let T: ¥ > X be a mapping on X and let X bea 
T-orbitally complete metric space. If there exists real numbers «;, 8 for every x, y © x 
such that «, + a, + a3 > 8 and — a, = 0 or B — as = O, and 


a, P [Tx, Ty] + %, P [x, Tx] + as P Ly, Ty] 
+ a, min {p [x, Ty], P [y, 7x]} < BP Lx, y] --(7) 
then for each x € Y, then sequence {7” x} converges to a fixed point & of T. 
Proof is analogous on the proof of the Corollary 6. 


Also, we note, Corollary 6 of Dhage* which the assumption of the continuity of 
T is removed. Also, the condition (6) implies out the condition (7). 


Let X denote a metric space and K a bounded subset of ¥. Following 
Kuratowski® we denote by « (K) the infimum of all « > 0 such that K admits a finite 


converging with subset of diameters less that e. We use the following properties of 
the number a. 


1, « (K) = Oif and only if K is precompact. For this reason « (K) is called the 
measure of noncompactness of K. 


2. « (AUB) = max {a (A), « (8)}, (A, B CX). 


3. If K is compact then «(K) = 0. Also, a (CIK) =O<=« (K) = 0, and 
0 S«(K) S38 (K). 


4. «(K) = 0 and ¥ complete imply K is compact. 


We list below some of contractive mappings for which various fixed point theorems 
have been stablished : 


(1) (Furi and Vignoli‘). The continuous mapping T is called densifying, if for 
every bounded subset A of X, such that x (4) > 0, we have « (T (A)) < « (A). Let 
F be a real lower semicontinuous function defined in ¥ x ¥. The densifying mapping 


Tis said to be weakly F-contractive if the condition F (Tx, Ty) < F(x, y) holds 
for allx,y € X¥,(x Hy). 
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(2) (Ray and Chatterjee", Singh’?, Khan’). Let F: ¥ x Y¥ > R{. be continuous 


and T : X + X be densifying mapping such that F (Tx, Ty) < aF (x, y) + BF (x, Tx) 
+ yF(y, Ty) for each pair, of distinct points x, »y € X and for nonnegative real 
numbers «, 8, y with « + 8 + y < 1, then T is called generalized densifying. 


(3) (Bose and Mukherjee") Let F be a continuous symmetric mapping of ¥ x X 
into the set of nonnegative reals, such that F (x, y) = Oiffx = y and 


F (Tx, Ty) < aF (x, Tx) F (y, Ty) (F (x, y))"! + bF (x,y), x fy 
for each pair of distinct points x, y € X¥ wherea + b < 1 (a,b > 0). 


Corollary 10 (Furi and Vignolit)—Let T: X > X bea densifying (therefore, by 
definition, continuous) and weakly F-contractive mapping defined in a complete 
metric space. If for some x, € X, the sequence {x,} of iterates starting from x, = x 
is bounded, then T has a unique fixed point & in X. 


Proor: For O(x)we have TO (x) C O(x), then O(x) is an invariant set. 
Denote by ClO (x) the closure of O (x). CIO (x) is an invariant set, too; indeed, from 
the continuity of T it follows T (Cl O (x)) C Cl T(O (x)) C Cl O (x). Now we shalt 
prove that ClO (x)is compact. For this it is sufficient to show that « (O (x)) = 0, 
since ina complete metric space the precompact sets are also relatively compact. 
Suppose « (O (x)) > 0; in this case 


a (T (O (x))) < « (O (x)) and since O (x) = TO (x) U {x} ...(8) 
it follows that 
a (O (x)) = max {« (TO (x)), « (x)} + max {« (TO (x)), 0} = «(TO (x)). 
But this contradicts (8), hence « (O (x)) = 0; 
Also, let B (x) = F(x, Tx). Since Tis continuous and Fis lower semiconti- 
nuous, B is lower semicontinuous, it is easy see to that B: ¥ > R}, satisfies all the re- 


quired hypothesis in localization monotone principle. Hence, applying the localization 
monotone principle we obtain T& = &, for some & € ¥. Uniqueness, it follows 
immediately from the weak F-contractivity. 


Analogous to the proof for the preceding statement we have the same proof for 
the statements of Ray and Chatterjee’’, Singh’*, Khan‘, Bose and Mukherjee’. 
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Recently the notion of feebly closed mappings was introduced. In this paper 
this concept is shown to coincide with the notion of «-closed mappings. 
Furthermore it is shown that if the codomain is appropriately retopologised 
the concept of feebly closed mappings coincides with the usual notion of clos- 
ed mappings. Some properties of «-closed mappings are investigated. 


1. INTRODUCTION 


Let S be a subset of topological space (¥, 7). We denote the closure of S and 
the interior of S with respect to F by Fel S and F int S respectively, although 
we may suppress the J when there is no possibility of confusion. We denote the topo- 


logy induced by F on S, by Fs. Njastad® introduced the concept of an «-set in 
(X, 7). Asubset S of (XY, F) is called an «-set if S C F(int Fel (F int S)). The 
notions of semi-open set and preopen set were introduced by Levine? and Mashhour 
et al.5 respectively. A subset S of (X,Z ) is called a semi-open set (respectively preopen 
set) if SC Gcel(Fint S) (respectively Sc Gint(Fel S)). The complements of 
an a-set, a semi-open set and a preopen set are called «-closed, semi-closed and pre- 


closed respectively. Following Nijastad® we denote the family of all «-sets in (X,F) 
Dy as Nijastad proved that F* is a topology on X. As any open set in (X¥, H) is 
an a-set, DC G* inthe lattice of topologies on the set X. If A isa subset of 
(X, 7), then the intersection of all semiclosed sets containing A is called the semi clo- 
sure of A, and is denoted s cl A. The largest semi-open set contained in A is denoted 
by s int A. Maheshwari and Tapi’ defined A to be a feebly open set in (X, 7) if there 
is an open set Usuch thatU C ACs cl U. The complement of a feebly open set is 
called a feebly closed set. 


2. Frepry CLosep MAPS 


The concepts of «-closed and feebly closed mappings have been introduced by 
Mashhour et al.§ and Maheshwari and Jain respectively. 


*The second author wishes to acknowledge the support of the University of Auckland Research 


Fund. 
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Definition 1—A functionf: (X, FZ) > (%, & ) is called 

(i) feebly closed if the image of each closed set in X is feebly closed in Y; 

(ii) «-closed if the image of each closed set in X is «-closed in Y. 

Lemma 1—Let A be a subset of (X, Z). Then sint (cl A) = cl (int (cl A)). 

Proor : Notice that cl (int (cl A)) is a semi-open set since cl (int (cl A)) = cl (int 
(cl (int (cl A)))), and cl (int (cl A)) C cl A. Therefore cl (int Cl A)) C s int (cl A). 


Conversely, if U is any semi-open set contained in cl A, then U C cl(int VU) C 
cl (int (cl A)) and therefore s int (cl A) C cl (int (cl A)). 


Maheshwari and Jain® in Lemma 3 of their paper showed that a subset A of 
(X, Z) is feebly closed if and only if s int (cl A) C A. 


Proposition |—If A is a subset of (X, 7), then A is feebly closed if and only if A 
is a-closed. 


Proor : It follows from the definitions of an a-set and an «-closed set that a 
subset A of (X, @) is «-closed if and only if cl (int (cl A)) C A. 
Since cl (int (cl A)) C A if and only if s int (cl A) C A, by Lemma 1, A is a-closed 
if and only if A is feebly closed. 

An alternative proof of this result has been given in Proposition | of Jankovic and 
Reilly’, who proved that if A is a subset of (Y¥, F), then A is feebly open if and only 
if A is an a-set. 


Proposition 2 follows immediately from Proposition 1 and Definition |. 

Proposition 2—The following are equivalent : 

(1) f: (X%,F)—> (Y¥, &) is feebly closed; 

(2) f: (X, FZ) > (¥%, &) is a-closed; 

(3) f: (X, Z) > (Y, Q°*) is closed. 

If the codomain space of a feebly closed mapping / is retopologised in an obvious 
way, then fis simply a closed mapping. This observation puts the notion of feebly 
closed mappings into a more natural setting, and enables us to provide immediate 


proofs of some of the results in Maheshwari and Jain®. For example, Propositions 5 


and 6 of Maheshwari and Jain® are well-known results (Murdeshwar*, Theorem 4.26, 
P. 96 and Theorem 4.28, p. 26) restated in this setting. 


3. a-CLOSED MAPPINGS 


The following classes of generalized closed mappings were introduced in 
Mashhour et al.° and Noiri?®. 
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Definition 2—A function f: (X¥, F) > (Y, W&) is called 

(i) semi-closed if the image of each closed set in X is semi-closed in ¥; 

(ii) preclosed if the image of each closed set in Y is preclosed in Y. 


It is shown in Theorem 3 of Reilly and Vamanamurthy” that a subset of (X, 7) 


is an a-set if and only if it is semi-open and preopen. Thus we have the following 
result. 


Proposition 3—A function f : (X¥, FH) > (Y, &W) is «-closed if and only if it is 
semi-closed and preclosed. 


Examples | and 2 show that the separate converses are not in general true. 


Example |—Let X = {a, b, c} and define the topologies FG to be the discrete 
topology and 27 = {¢, x, {a}, {c}, {a, c}}. We define f: (X, FZ) > (X, WY) by f(@ 
= f(b) =f(c) =a. Then fis preclosed but not «-closed since {a} is preclosed in 
(X, 27) but not «-closed in (X, 97). 


Example 2—If Fis the discrete topology and 9 is the indiscrete topology then 
tf: (% FZ) —- (% &) is semi-closed but not «-closed. 

Andrijevic! showed that if M is a subset of (¥, Z) then (Fm)* C (F*)m (his 
Theorem 3.2) and if M is preopen then (F m)* = (F*)m (his Theorem 3.6). 

Proposition 4—Let f: (X¥, FZ) > (Y, WM) and f(X) CY: C Y. 

(1) If f: (X%, F) > (%, QW) is a-closed, then so is f:(%, F)>(% &). 


(2) If Y, is preopen and f: (X, W) — (Y,Q7) is «-closed, then f: (X,F) > (%, 
Wy) is z-closed. 


Proor : (1) If f: (X, FZ) > (%, Wy.) is -closed then f: (X, F) > (Vi, (Wy, )*) 


is closed. By Andrijevic’s result f: (X, FZ) > (Y%,, (YH *)y,) is closed and by Theorem 


4.24 (2) of Murdeshwar' f:: (Y, FZ) — (Y, Q*%) is closed. Hence f: (X, FZ) > (YQ) 
is a-closed. 

(2) Iff:(X%, TZ) > (Y, Q) is a-closed then f : (X, TF) (Y, QW") is closed 
and by Theorem 4.24, of Murdeshwar,’ f:(X, 7) > (hh, (2*)y,) is closed. 


If Y; is preopen then by Andrijevics result f:(X, F) > (%1, (W1)*) is closed, and 
sof: (X,F7)7(%, Wy ,)is a-closed. 


Proposition 5—If f is an o-closed mapping f: (¥, 7) > (¥, W@), A: C F and 
* CW, then f: (X, Fi) > (Y, W) is a-closed. 
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Proor: Iff: (1. F)> &% Dis a-closed then f:(X, 7) > (Y¥, @™*) is 
closed, and by Theorem 4.23 of Murdeshwar*, f: (X, Fi) > (Y, W; ) is closed. Hence 
fel Xe sy) ee Ul , W,) is «closed. 

Proposition 6—If f : (X, FZ) > (Y, W) is an a-closed mapping and B, G Cc ig 
then if f-! (B) and f-! (C) have disjoint neighbourhoods, B and C have disjoint neigh- 
bourhoods in (Y, Q«). 

Proor: Iff:(X,F) > (%,Q) is a-closed then f: (X, F) > (Y, YH") is 
closed and the result follows from Murdeshwar? Theorem 4.28 (2). 


Proposition 7—Let f : X > (Y, Q@) and let X be given the preimage topology 7. 
Then / is a-closed if and only if f (xX) is «-closed in (Y, Q/). 


Proor : One implication, namely if f is «-closed then f (X) is «-closed in (Y, 27) 
is clear since Y is closed in the preimage topology. 


Conversely, if f (X) is «-closed in (Y, Q/) thenf(X) is closed in (Y, Q/ Sh Lat 
F,, be a preimage topology on X for f: X¥ > (¥, W*). 


Then, by Theorem 4.30 of Murdeshwar' f:(X, Fa) > (Y, Q/*)is closed. Since 
W C W* and therefore F C Ga, f: (X, FT) > (Y, W*) is closed, and so f: (X,F) 
— (Y, Q&) is «-closed. 


It is well known that the 7; property is preserved under closed mappings. The 
following example shows that the 7, property is} not preserved under «-closed 
mappings. 


Example 3—Let X be an infinite set and p be a fixed point of ¥. We definea 
topology Gon X as follows: forG CX,GE FS ifG=¢or G=XorxX —G is 
finite. We define a topology Qf on X as follows: for each GCY,GEQO” if 
(i) G= ¢or G = X, or (ii) G C X¥ — {p} and X — Gis finite. (Y, F) is T, but (¥, @) 
is not T since for any point x distinct from p the only open set containing p, namely 
X, contains x. Let f : (X,F) —(X,Q/) be the identity function. Then / is «-closed since 
if A is a closed subset of (X, F) then either A is a closed subset of (X, 27) and there- 
fore «-closed in (X, Q/), or A is finite, nonempty and p & A. In this case O cl (@ int 
(@ cl A)) = ¢ C A and therefore A is an «-closed subset of (Y, QW). 


The following proposition is a generalization of the well known result that regu- 


larity is preserved under continuous, open and closed surjections (Murdeshwar* Theorem 
12.14, p. 206). 


Lemma 2—If U and V are subsets of (X, 7), U is open and UC V, then 
cl U C cl (int (cl V)). 


ProoF: If U is open and U C V, then U C cl V so that U C int (cl V). There- 
fore cl U C cl (int (cl (V)). 
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. Proposition 8—If f: (X, FZ) > (Y, Q/) is an open, continuous, «-closed surjec- 
tion and (X, 7) is regular, then (Y, Q/) is regular. 


Proor : Let p € Y and U be an open set in (Y, 97) containing p. Letx € X 
such that f(x) = p. Since (Y, @) is regular there is an open set V in (X, FZ) such 
thatxE VC GealVCf'(U) sothat pEf(VV) Cf(FealV) CU. Since fis «- 
closed, f(F cl V) is «-closed and since f is open f (V) is open so that, by Lemma 2, 
LAfV) C Wel (& int (WM cl f(Hel f (V))) CU and therefore p € f(V) 
C Hel f(V) C U. 

It is well known that normality is preserved under closed, continuous surjections 
Murdeshwar* Theorem 15.3 (1). The following proposition is a generalization of this 
result. 


Lemma 3—If U and V are subsets of (X¥, 7) and UM V = ¢ then int (cl (int 
U)) CM int (cl (int V) = ¢. 

Proor : If U(\ V = ¢thenint U int V = ¢, so that int UN cl (int V) = ¢. 
Therefore int U (-) int (cl (int V)) = 4 which implies that cl (int VU) int (cl (int V)) 
¢, so that we have int (cl (int U)) CM int (cl int V)) = ¢. 

Proposition 9—\f f:(X, FZ) > (Y, W) is a continuous, «-closed surjection and 
(X, 7) is normal, then (Y, 27) is normal. 

Proor : Let A and B be closed sets of (Y, 27). Then there are open disjoint 
sets U and V in(X, G) such that f-!(A) C U and f7*(B) CV, by normality of 
(X, F). By Proposition 6 there are disjoint «-setsC and Din (Y,Q/) such that 
AC Cand BC D>s0 that, by Lemma 3, int (cl (int C)) and int (cl (int D)) are dis- 
joint open sets in (Y, 2/) containing A and B respectively. 
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In this paper, attempts have been made to generalize the concept of set- 
connected mappings in bitopological spaces. Such a mapping has been 
characterized and some of its properties have been studied. Finally, its 
relationships with pairwise continuous and pairwise weakly continuous 
mappings have been investigated. 


1. INTRODUCTION 


Kwak? introduced ina topological spacea new type of mapping, called set- 
connected mapping which is weaker than a continuous map. Noiri® further investi- 
gated different properties of this mapping and showed that every weakly continuous 
surjection is set connected, but not conversely. In this note, we have studied the 
corresponding results in the bitopological setting. Since a bitopological space, as 
introduced by Kelly’ as a natural structure, is a richer structure than that of a topo- 
logical space, it is of much use and importance to study the generalizations of topo- 
logical notions and implications in bitopological situation. Apart from this, it is 
interesting to notice the interactions of the two topologies in a bitopological space 
with reference to a study of a bitopological notion. These are the underlying reasons 
for extending and studying set connected maps in a bitopological space. 


Section 2 opens up with the study of different properties of p-connected mapping 
(henceforth, the word “pairwise” shall be abbreviated as “‘p—’’), specially in relation to 
its inherent applications to p-connected spaces, as introduced by Pervin’. Relations of 


this mapping with p-continuous and p-weakly continuous mappings have been discussed 
in section 3 and 4 respectively. 


Throughout this note Y and Y always represent bitopological spaces (XY, Q,, Q,) 
and (Y, P,, P,) respectively and f: X — Y denotes a mapping from X into Y. If ACX, 
then the subspace A is the bitopological space (A, (Q,)4,(Q,)4) induced by 
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(X, Q:, Q.). Q:—cl A and Q)-int A respectively shall denote the closure and interior of 
A in (X, Q:), fori = 1,2. An ij-coset in any space (X, Q,, Q2) is a Qi-closed and Q,- 
open set in X and i, j can assume values | and 2 with i + /. 


2. P-sET CONNECTED MAPPINGS 


Definition 2.1—A bitopological space ¥ is said to be 12-connected (21-connected) 
between A and B (where A, B are non-empty subsets of X) iff there exists no 12-coset 
(2l-coset) F such that dA C FC X — B. 


Clearly, X is 12-connected between A and B iff it is 21-connected between 
Band A. 


Definition 2.2A—X is said to be p-connected between A and B (where A, B are 
nonempty subsets of X) iff X is 12-connected as well as 21-connected between 
A and B. 


Remark 2.3: If X is ij-connected between A and Band C > A, DD B, then X 
is ij-connected between C and D. 


Definition 2.4—A mapping f: X > Yis said to be p-set connected iff f(X) is 
ij-connected between f(4) and f(8) in the subspace bitopology whenever YX is ij- 
connected between A and B, where i, / = 1, 2 andi + j. 


Remark 2.5: It is clear that the definition 2.4 can equivalently be modified by 
writting 12 or 21 in place of ij. 


Lemma 2.6—If a subspace M of a bitopological space X is /j-connected between 
A and B then so is the whole space (where, as usual, i, j = 1, 2 andi # /). 


Proor : If there is any ij-coset F in the whole space X with A C FC X—B, 
then F (-) M is ij-coset in M with AC FQ M CM — B. 


Lemma 2.7—A subspace M which is 12- (or 21-) coset in X, is ij-connected 
between A and B if X is ij-connected between two subsets 4 and B of M. 


Proor : If Fis any 12-coset in M, then F is also 12-coset in X, as M is 12-coset 
in X. This fact implies that Mis 12-connected between A and B if X isso. Now 
if X is 21-connected between A and B then it is 12-connected between Band A; so M 
is 12-connected between B and A which implies that M is 21-connected between A and 
B. Similar argument is applicable when M is 21-coset in X. Hence the lemma is 
proved. 

Theorem 2.8—A mapping f: X > Y is p-set connected iff f-! (F) is 12-coset in 
X for any 12-coset F in f (X). 

Proor : Let f be p-set connected and F be any 12-coset in f (X). If f~? (F) is not 
12-coset in ¥ then X is 12-connected between f~' (F ) and X -f-!(F) and therefore, 
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as fis p-set connected, f (X) is 12-connected between f(/ -1(F)) and f(X — f7™ (F)). 
But f(f> (F) = FOS (X) = F and f(X — £7 (F)) = f(X) — F imply that Fis 
not 12-coset in f(X) — a contradiction. Hence f~! (F) is 12-coset in X. 


Conversely, let the condition hold for fand let X be 12-connected between A and 
B. If f(X) is not 12-connected between f (A) and f (B), then there exists a 12-coset 
Fin f(X) such that f(A) C FC f (X) —/S (B). But we find that A Cult CA) 
Cf (F), f-1(F) C X — B and f~) (F) is 21-cosetin X. This contradicts that X is 
12-connected between f (A) and f (B) i.e., f is p-set connected. 


Remark 2.9 : (a) Similarly as in the above theorem, we can show that f is p-set 
connected iff f-! (F) is 2l-coset in X for any 2l-coset Fin f(X). (b) If/f is p-set 
connected, then f-! (F) is 12-coset (21-coset) in X for any 12-coset (21-coset) F in Y. 


Definition 2.10’—A bitopological space (X, Q,, Q.) is called pairwise connected 
(p-connected, according to our abbreviation) iff ¥ cannot be expressed as the union of 
two non-empty disjoint sets 4 and B such that (A ™ Q,—cl B) U (Q,—cl A MB) = ¢. 


It is shown’ that (XY, Q, Q,) is pairwise connected iff X¥ cannot be expressed as 
the union of two non-empty disjoint sets A and B such that A is Q,-open and B is 
Q.-open. 


Now Lemmas 2.11 and 2.12, stated below, follow immediately. 


Lemma 2.11—Every mapping f: Y¥ > Y such that f (X) is p-connected, is p-set 
connected. 


Lemma 2.12—Let f : X > Y be a p-set connected mapping. If X is p-connected, 
then f (XY) is p-connected. 


. Lemma 2.13—Let f: X > Y be p-set connected and A (C X) be such that f (A) 
is 12-coset (or 21-coset) in f(X). Then the restriction f|A: AY is p-set connected. 


Proor: Let A be 12-connected between C and D. By Lemma 2.6, X is 
12-connected between C and D, so f (X) is 12-connected between f (C) and f(D). As 
f(A) is 12-coset (or 21-coset) in f (X), by Lemma 2.7, f(A) is 12-connected between 
f(C) and f(D). Hence by Remark 2.5, the lemma is proved. 


Theorem 2.14—Let f:X—Y be p-set connected, by-open (i.e., £:(X, Oy) > 
(Y, Pi) is open for i = 1, 2), surjection and f~* (y¥) be p-connected for each y of Y. 
Then for any 12- (or 21-) coset Fin Y, F is p-connected iff f-1 (F) is p-connected. 


: Proor : We prove the theorem for the case when Fis 12-coset in ¥. Similar 
will be the proof when F is 21-coset. Let f~ (F) be not p-connected for some 12-coset 
Finy. Then there exists a 12-coset K with ¢# KAf—(F), in the subspace 
bitopological space (f~* (F), Q,/f-} (F), Q,/f-1 (F)). We show that f (K) is 12-coset 
in F with P$AL(K)JAF. Asf(y) is p-connected, either f-} OC K, or f~ (y) 
Cf (F) — K, for all y € Fand so f(K) # F (+ 4) and f (K) Afgn (F) — K) 
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=. As f is surjective, f(K) U f(f-'(F) — K) =F and so f(f- (F) — K) 
= F — f (K). Now i being bi-open f/f~!(F) is bi-open onto F and hence f (K) is 
12-cosetin F. This implies that F is not p-connected. So Fis p-connected implies 
f- (F) is p-connected. 


Next, since f (f~! (F)) = F which is 12-coset in Y, by Lemma 2.13, the restric- 
tion f/f-* (F): f~ (F) > Y is p-set connected. Now by Lemma 2.12 it follows that 
Sif (F) (f7 (F)) = Fis p-connected if f~! (F) is p-connected. Hence the theorem. 

Theorem 2.15—Let f: X > Y be a p-set connected mapping. Then for each 
P € f(X), f(Q: — cl f-!(p)) C Cp, where C;, is any ij-coset in f(X) (or in Y) contain- 
ing the point p (where, as usual, i, 7 = 1,2 andi ~/). 

. Proor: As C, is ij-coset in f(X) (or in Y) and f is p-set connected f~ (C,) is 
ij-coset in X (by Theorem 2.8 and Remark 2.9 (b)). Also f-!(p) C f7! (C,), therefore 
Q:—cl (f"(p)) C f7(C,). Hence f(Q: — cl (f-*(p)) C ff-* (Cr) C Cp and 


the theorem is proved. 


Theorem 2.16—Let (X, Q,, Q2) be a bitopological space. If a sequence {p,} 
converges to p in (X, Q;) and also if phas a p-connected Q)-open neighbourhood, 
then there exists some , such that Yis p-connected between p and p,, for all n > 1% 
(i =f, 2). 

Proor: If pand p, are contained ina p-connected set then X is p-connected 
between p and p,. 


Definition 2.17'—Let (X, Q,, Q.) be a bitopological space. Q: is said to be 
locally connected with respect to Q, if for each point x in X and each Q,-open neigh- 
bourhood U of x there is a pairwise connected Q,-open set G such that x € GC U. 
(X, Q:, Q;) is pairwise locally connected if Q, is locally connected w. r. t. Q, and Q, is 
locally connected w.r.t. Qi. 


Corollary 2.18—If a sequence {x,} converges to x in (X, Q1) and Q; is locally 
connected w.r.t. C; then there exists some Mo such that ¥ is p-connected between 
p and p,, for all n > ng, where i = 1 and j = 2, ori = 2 andj = 1. 


Definition 2.193—In a bitopological space (X, Q,, Q.), Qi is said to be extremally 
disconnected w.r.t. Q; iff Q;-closure of any Qi-open set is Qi-open (i, j = 1, 2 and 
i+j). The space is p-extremally disconnected iff QO, is extremally disconnected 
w.r.t. QO. and Q, is extremally disconnected q.r.t. Q1. 

Definition 2.20‘'—A_bitopological space (X, Q:, Q2) is said to be pairwise 
Hausdorff if given distinct points x, y of X, there is a Q,-open set U and a Q2-open 
set V, such thatx € U,yEV,UN V=4¢. 

Theorem 2.21—Let f : X + Y be p-set connected and Y be p-Hausdorff and P; 
is extremally disconnected w.r.t. P; in Y. Then the graph G(/) of fis closed in 
(Xe 7, 0, xX Pi) GI = 1,2 andi +‘). 
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Proor : Let (x, y) & G(/) so that y A f(x). Then y has a P;-open neighbour- 
hood whose P,-closure, say V, does not contain f (x) (since Y is p-Hausdorff). As P; 
is extremally disconnected w.r.t. P;,V is /i-coset in Y and so f~ (V) is ji-coset in X 
[by Remark 2.9 (b)] and x & f7'(V). LetU=X —f4(V). Then U is Q,-open 
in X containing x and f(U) 1 V=¢. Thus U x Vis an open neighbourhood of 
(x, y) in (X x Y, Q; x Pi) with (U x VIN G(S) = ¢. 
Hence G (/) is closed in (¥ x Y, Q; X P;). 


Remark 2.22: In the above theorem, one might expect that G (/) is closed in 
(X x Y, O; x Pi) (fori = J, 2). That this is false is shown by the next example. 


Example 2.23—Let X¥ = Y = R = real line. Let Q: = P,; = cofinite topology 
on R and Q, = P, = discrete topology. Letf: X > Y be the identity map. Then 
f is clearly p-set connected, (Y, P, P,) is p-Hausdorff and P, is extremally disconnected 
w.r.t. P,. But G(f) is not closed in (X x Y,Q, x P,). In fact, G (f) is dense in 
(x xk Oi exatak 


Theorem 2.24—Let f: X — Y be a mapping andg: X — (X x Y,Q, xP,, Q, 
x P2) be given by g (x) = (x, f(x)), forx € ¥. Then fis p-set connected if g is so. 


Proor : Let F be 12-coset in f (X) C Y. Then X x F is 12-coset in the subspace 
X x f(X) C X x Y(we note that the topologies on X¥ x f(X), namely Q;x Pi/ f (X) 
and that inherited asa subspace from Q; x P; are identical, for i= 1,2). As 
g(X) C X x f (X), (X X F) M g (X) is 12-coset in g (X) andso g7' [(X x F) ) g(X)] 
is 12-coset in X by Theorem 2.8. But g"[(X¥ x F)( g(X))}=27(X x F) =f- 
(F). Hence / is p-set connected. 


Remark 2.25 : Converse of the above theorem is false as is seen from example 
below. 


Example 2.26—Let X = Y = R, the real line; Q, = P, = the left hand topology 
and Q, = P, = the right hand topology. Let f: X — Y be defined as S(x)=~ if x30 
and f/(0) = — 1. Here fis p-set connected as Y has no nontrivial 12- or 21-coset. 
Now [0, oo) x [0, 00) is closed in Q, x Pi, (O,c0) x (0,e°) is open in Q, x P, and ([0,ee)) 
x [0,e°)) 1 g (X) = (0,00) x (0,00)) A g(X) = K (say). Thus K is 12-coset in the 
subspace g (X) of X x Y. ButX has no non-trivial 12-coset whereas g 1 (K) = (0, oo). 
Hence g is not p-set connected. 


3. P-seET CONNECTED MAPPINGS AND P-CONTINUOUS MAPPINGS 


. Definition 3.1’—A mapping f: X + Y will be said to be p-continuous iff the 
induced maps f : (X, Q,) > (Y P,) and f : (X, Q,) > (Y, P.) are continuous. 


By virtue of Theorem 2.8 it follows that : 


Theorem 3.2—Every p-continuous mapping is p-set connected. The following 
example shows that the converse is false. 
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Example 3.3—Let X = Y = R, the real line, and P; = Q, = the usual topology 
on R and P, = Q, = the countable complement topology on R. Let f: X > Y be 
defined by f(x) = x, x ~ Oand f(0) = 1. Then /is p-set connected, since there is 
no proper 12- or 2l-coset in f(X) = Y — {0}. But f:(X,Q,) > (Y P,) is not 
continuous at x = 0 and hence f: X¥ — Y is not p-continuous. 


Now, we proceed to find the conditions under which a p-set connected mapping 
is p-continuous. 


Theorem 3.4—Let f: X + Y be p-set connected. If Y be p-Hausdorff and P, 
(or P,) is extremally disconnected w.r.t. P, (or P,), then f/C : C > Y is constant for 
every p-connected subset C of X. 


Proor: Let x,y & C andx+~y. If possible, let f(x) 4 f(y) in Y. Then 
there exists 21-coset V in Y such that f(x) € V and f(y) & V, since Y is p-Hausdorff 
and P, is extremaly disconnected w.r.t. Ps. Now, f~'(V) is 21-coset in X as f is p-set 
connected and therefore, f-! (V) M C is a non-empty proper 21-coset in the subspace 
C. This condradicts that C is p-connected. Hence f(x) = f(y), for allx,y EC, 
i.e., f/C : C — Y is constant. 


Corollary 3.5—Let f: X > Y be p-set connected, Y be p-Hausdorff and P, 
(or Ps) is extremally disconnected w.r.t. P, (or P:) and for each x € X it has a 
Q,-open p-connected neighbourhood. Then f~? ( y)is Qi-open in X for each y € ¥, 
te a) Ee = 


Proor: Follows from the fact that for each x € f-'(y) there is a Q;-open 
p-connected set U with x € UC f? (7): 


Remark 3.6: In the above corollary, if for each x € X, x has a Q;-open 
p-connected neighbourhood and a Q,-open p-connected neighbourhood, then f is 
p-continuous Thus 


Corollary 3.7—If f: X — Y be p-set connected, XY be p-locally connected and Y 
be p-extremally disconnected and p-Hausdorff, then / is p-continuous. 


Definition 3.8—A space (X, Q,, Q,) is p-c-compact iff every cover of any QO: 
closed set in ¥ by Q,-open sets of X has finite subfamily whose union is Qj-dense in 
the Q;-closed set, for i, 7 = 1,2 andi + /. 


Theorem 3.9—Let Y be p-extremally disconnected, p-c-compact and p-Hausdorff. 
Then f: X — Y is p-continuous if it is p-set connected. 


Proor: Let fbe not p-continuous. Then there exists a Pj-closed set F in Y 
such that f~ (F) is not Q;-closed in X (fori = lor 2). Letx € Qi-el(f” (F)) 
—f-(F). Then X is ij-connected between f~! (F) and x. So f (X) is ij-connected 
between f (f~! (F)) and f(x). By Lemma 7.6 and Remark 2.3, we obtain that Y 1s 
ij-connected between F(> ff" (F)) and f (x). Now as Y is p-Hausdorff, for each 
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y € F there exists a P;-open neighbourhood Vy of y in Y such such that f (x) & P,—cl 
(Vy). Then the family {V): y € F} is a cover of F by P,-open sets inY. By 
p-c-compactness, there exists finite set of points V1, Y2,---,¥n In F such that F C Pi—cl 


( U V,) = V(say). Then V is ij-coset in Y, since Y is p-extremally disconnected and 
r= ’ 


Pi—cl((U Vy) =U Piel (Vy,): Also, f (x) & V, since f (x) & Pi-cl (Vy), for any 
dal | r rol 


y € F. This contradicts that Yis ij-connected between Fand f(x). Hence f is 
p-continuous. 


4. P-seT CONNECTED MAPPING AND 
P-WEAKLY CONTINUOUS 
Definition 4.1°—A mapping f : (X, Q:, Q.) — (Y, P: P.) is said to be Q, Pi-weakly 
continuous w.r.t. P, if for each point x of X and each P,-neighbourhood V of f (x) 
there exists a Q,-neighbourhood U of x such that f(U) C P,—cl (V). 


The mapping / will be called p-weakly continuous if it is Q, P: weakly continuous 
w.r.t. P, as well as Q, P, weakly continuous w.r.t. P, 


Lemma 4.2°—A mapping f: X > Y is Q, P,-weakly continuous w.r.t. P, iff 
for each P,-open set V of Y, f-? (V) C Q,—int (f= (P,—cl V)). 


Lemma 4.3?—If a mapping f: X > Y is Q, P; weakly continuous w.r.t. P, then 
Q,--cl (f7' (V)) C f-! (Pi—cl V), for any P,—open set V. 


Theorem 4.4—If a surjection f : (X, Q,, Q.) > (Y, Pi, P2) is p-weakly continuous, 
then f is p-set connected. 


Proor: Let V be 12-coset in ¥(= f(X)). Since V is P,-open, by Lemma 4.2 
f~ (VV) C Q.—int (f-1 (P,—cl V)) = Q.—int (f -' (V)), since V is P,-closed. This im- 
plies f~' (V) = Q, ~ int (f~' (V)) and hence f~! (V) is Qo-open in X. Again by Lemma 
4.3, 0,—cl(f (V)) C f+ (P,-cl V) = f-1 (V) so that 0:-—cl(f-!(V)) =f (Vv) 
and hence f~! (V) is O.-closed in X¥. Thus f~ (V) is 12-coset in ¥. Then by Theorem 
2.8, fis p-set connected. . 


Remark 4.5 : The following example shows that in the above theorem the condi- 
tion “‘p-weakly continuous” can not be relaxed by Q, P,-weakly continuous w.r.t. P, 
or Q, P.-weakly continuous w.r.t. P;. ; 


Example 4.6—Let X = Y = R, the real line, Q1, Q, be respectively the discrete 
and indiscrete topology on R; P, the lower limit topology and P, the usual topology 
on Y. Suppose f: (X, Q,, Q.) > (Y, P,P.) be the identity mapping. Now A = [0, oo) 
is 21-coset in Y. But f-!(4) = Ais not 21-coset in X. Therefore, Jf is not p-set 


connected. But itis clear that fis Q, P; weakly continuous w.r.t. P, and is not 
p-weakly continuous. > 
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Remark 4.7: The converse of the Theorem 4.4 is false as is seen from the 
example below. 


Example 4.8—Let ¥ = Y = R; Q,, Q, respectively denote the cofinite and 
discrete topology on X and let P; = P; = usual topology on Y. Let f: (X, Q:, Qs) 
— (Y, Pi, P2) be the identity map. Since (Y, P:, Pz) is p-connected, f is p-set connected. 
Now consider o € X¥ and V = (— I, !)isa P,-neighbourhood of f{o) = 0 in y; 
Then there cannot exist any Q,-open neighbourhood U of o such that f(U) C P,.—clV 
=[-—1, 1]. Hence fis not p-weakly continuous. 


Theorem 4.9—Let f: X > Y be p-set connected and Y is p-extremally disconne- 
cted, then f is p-weakly continuous. 


Proor: Let x € Xand f(x) EVE P, (= 1,2). Then P,—clV is ji-coset in 
Yandso (P;—clV)  f(X) is ji-coset in PAO.OT CRN ge PRP a a Since / is p-set 
connected f-' (Pj;—cl V 1 f (X)) = U (say) is ji-coset in Y. Then U is a Q)-open 
neighbourhood of x such that f(U) C P;—cl V (i,j = 1, 2 andi # j), Hence fis 
p-weakly continuous. We shall sharpen the result in Theorem 4.9. For that we 
require the following definition. 


Definition 4.10°—A mapping f : (X, Q1, Q.) > (Y, P,, P2) is called Q, P,-almost 
continuous w.r.t. P, iff for each x of X and each P,-neighbourhood V of St (x) in Y 
there exists a Q.-neighbourhood U of x in X such that f(U) C P,- int (P,—cl V); «If f 
is, in addition, Q, P,-almost continuous w.r.t. P,, then f is called p-almost continuous. 


Remark 4.11: Since every p-almost continuous map is obviously p-weakly 
continuous, from Theorem 4.4 it follows that a p-almost continuous surjection is p-set 
connected. 


Now we observe that p-weakly continuity in Theorem 4.9 can be replaced by 
p-almost continuity. 


Theorem 4.12—Let Y be p-extremally disconnected and f: X > Y be p-set 
connected. Then/ is p-almost continuous. 


Proor : The same proof of Theorem 4.9 can be carried over here, noticing that 
P,—cl V=P;—int (P,—cl V) i, j=1, 2 and i $ /), since Y is p-extremally disconnected. 
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We define the class K, («; 4)—the class of «-starlike functions with respect 
to m symmetric points— consisting of f (z) = z + az” + .--- satisfying 


LOL 29 in pana SEO) +0 = 2f' 
azf’ (z)+ (1 — 2) fa (2) 


where A(z) is a given convex univalent function in E with / (0) = 1 Re A (z) 


<h(z), «70, 


> Oin E and fy (2) = "5" 0-4 f (wz) with @ = exp Qnijn), In this 
j=0 


paper we prove theorems which establish that K,, (a; h) C Kn (0; 4), Kn («;h) 


is closed under an integral operator and the coefficient estimate for this class. 


Also we define two class Cy, (x; 4)! of functions and investigate the proper- 
ties of this class. 


INTRODUCTION 


Let E = {z:|z|< 1} be the open unit disc in C, H(E) denote the class of 
functions f(z) holomorphic in E. Let A = {f € H(E):f(0) =f’ (0) — 1 = 0}. 
For a given positive integer n let w = exp (27i/n) and 


n— 


f. (2) = ee wo-I f(z), z © E. (1) 


Here we introduce new subclasses of S-the class of normalized univalent functions 
and study certain properties of these classes. First let us define the class K,, (a; A) 
which unifies the classes of convex univalent functions with respect to symmetric 
points and starlike univalent functions with respect to symmetric points. 


Definition 1--The function f(z) € A and bs st # Oin Eis said to be 


a-starlike with respect to n-symmetric points if it satisfies 
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az (zf’ z))’ + (1 — «) zf’ (z) 
az f’ (2) + (1 — alf,(z) 


=h(z), fore >0 


where A (z) isa given convex univalent function in E with h (0) = 1, Reh(z) > 0 
and ~ means subordination. Let us denote the class of such functions as K,, («;h). 


r 5 this reduces to the class K(x) of «starlike 


functions studied by Pascu and Podaru’. Also K, (1, i ) coincides with the 


When n = } and h (z) = 


(fara 
class Cs-defined by Das and Singh', K, ( 0, aes) ) is the class of starlike functions 


with respect to symmetric points investigated by Sakaguchi’. Of course K, (0;/) is 
the class of starlike functions with respect to n-symmetric points of Mocanu®. Also here 
we define another class C, (a; h) of x-close-to-convex functions with respect to n-sym- 
metric points. 

Definition 2—A function f(z) € A such that BOMAGK 0 in E is said to be 
a-close-to-convex with respect to n-symmetric points if it satisfies 

az (zf' (z))’ + (1 — «) zf’ (z) 
azg, (z) + (1 — 2) dn (2) 


= h(z), fora > 0 


a 


Ls 


where ¢, (Zz) = = 


Sw! ¢ (w'z) with 4 (z) © K, (a; h) andh(z) is as defined in 
j=0 


Definition 1. 


This class generalizes the classes defined by Pascu® and Das and Singh?. Let us 
investigate these two classes in this paper. In establishing our main results regarding 
these classes we often make use of the following lemma due to Padmanabhan and 
Parvatham‘. 


Lemma A—Let 8B, y € C,4 € H(E) beconvex univalent in E with h(0) = 1 
and Re (8h (z) + y) > 0, z € Eand let g € H(E) with qg (0) = 1 andg(z) x h(z), 
zEeE. If p(z)=1+p),2z +... in analytic in £, then 


p@) + gE) -« nie) > p@) 2h©) 





Now, we prove a lemma which we use in the sequel. 


Lemma—Let f(z) € K, («; A). Then f, (z) defined by (1) is in Ky («; 4). Further 


2f, (2) 


f, (Zz) 


zh (z) in E, 
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Proor : First let us consider 


az(zf, (z))’ +0 — alone (z) = \. "S az (zf (w! z))’ + (1-2) zf’ (w/z) 
az f(z) + (I — 4) fy (2) az fa(z) + (I— 2) fn (2) 


AS) az (2f' (!z)Y + (1 — a 2zf’(w%z) 
jae azf (z)+(l—« f, (2) 


Since f(z) € K, (x; 4) each term under the summation on the right side of the above 


equality is subordinate to h (z). Hence there exist 0, s in E such that 


oz ( af, (2) ) +(1'— a) zf- (2) 


=1 3 hg) =k) 
azf, (2) + (1 — af, @) ine 


> € E since h(E) is convex. Thus /, (z) € Ki (a, A). 
Since, 
az(of, @) +0- 2) of, @ 


a eas ee tae 2 p y+ azp’ (z) 
azf' (z)+ (l—a)f, (2) G (L'a) + api(z) 


where p (z) = of, (2) wee 
F.@” an application of a result by Eenigenburg et al. now yields 


af’, (2) 
ia p(z) = h(z). 
Theorem 1—We have the following inclusion relation : 


K, («; h) C K,, (0; A). 





e “ : tn ) 
PRrooF : Let f(z) € K, (a; h), p(z) = ao and qg(z) = : 


Then 
az (2f" (z))' + (1 — a) of’ (2) 
azf, (+ — fF, (2) 


(equation continued on p. 1117) 
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az(f, @PAa+/ALA+d -—wAP@sA 


= 


arh. WedetGl = 2) J, .(2) 
eefis (z) 
ce azp’ (z) + Fay, p{z) +(1 — «)p(z) 


Faiz) 
Fi ae 





aZ 


azf! (2) 
op” (2) +p) la -) + 
xz f' (2) 
(1-9) + 
azp' (z) 


Pita) baat Ge + p(z) 2 h (z) 


since f(z) € K,(«;h). Here q (z) ~ A(z) by the above lemma. Now an applica- 
tion of lemma A gives p(z) ~ / (z) in E which completes the proof of the theorem. 





—2z 
we geta result of Pascu® as a 


Remark 1: Whenn = | andh (z) = = 


particular case of our result. 


Now, we establish a theorem which shows that this class is closed under an in- 
tegral operator. 


Theorem 2—F (z) = iea-| rIe)-2 f(t) dt € K, (a; A) 
0 


whenever f(z) € K, («;h) for « > 0. 
Proor : Consider F, (z) defined by (1) with F (z) in the place of f(z). That is, 


n—1 ] z : 2! 
F,(z) = >. Zz w7! F (w/z). It is easy to see that F, (z) = spare J HE =* F{t)at. 


Differentiating this with respect to z and simplifying we get 


azF’ (z) + (1 — a) F, (z) = f, (2). te) 
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From the definition of F(z) we get, «zF’ (z) + (1 — «) F (z) = f(z) which on dif- 
ferentiation gives 
az. (ZF! (z))*-+ (Kia) 28’ (2) af (2): At), 
From (2) and (3) we get 
az (zF’ (z))’ + (1 — «) zF’ (z) rT 4-8 2 HD) 
azF’ (z) + (1 — «) F,(2) In @ 





because of Theorem 1 and the fact that f(z) € K,(a;h). This implies that 
F (z) € K, («;h). 


We now obtain the coefficient estimates for the class K, («;h). 


oo 
Theorem 3—Let f(z) = z + 2 a z'€ K, (a;h). Then ¥k > 1,n = 12 


tA tal | Ay | 


...(4) 
for i satisfying (kK — 1) n+1<i< kn+ L, WR alone Loe 


i (ia + 1 — a) 


where h(z) = 1 + h,z + .... 


; az (zf’ (z))’ +(1 — 2) 2f’ 
Proor :Let—— 2)" UL = &) 2” (2) =P(Z)=l+pz+p,224... 


azft, (z) + (1 — «)f, (2) 


l a=] 


where f, (z) = - 3 w-! f (wi z) and so fi) 2 - = f' (wiz). 


nN j-9 


as 2 L N ; 


az (zf" (z))’ + (1 — a) zf’ ()=(1+ p,z + P2z? + ...) (azf* (z) 


(l=) a) f(z) 


(equation continued on P. 1119) 
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=(l+p,z+...)(z +3q(¢#¢—0 4" ) > w(t=D)9) 71), 
in s=0 


Comparing the coefficients on either side, 


ja, | Ga + il — 2) S| A] + | Maly | (44 ) Sot 








+ | Alas | (=> ") ooeatt 


j=0 


eee thet) 


n 


n—-1 
x EZ w(t-2 
j—0 


+ | a | (Gert) "S w(t -2)F 


j=0 


(5) 


n—1l 
Also we make use of the fact that 2 of-DI = n when i= kn +1, ¥k>1; 
py 


n-l 
= wt) = 0 when isx~kn+ 1, ¥kS>1. 


470 
Hence fet us consider the two cases separately. 


Case (i) —First let us consider the casei = kn + |. Then we establish (4) for 
all k > 1, by the method of induction, In this case (5) reduces to 


kn (kna + 1) | Gknai | < | Ay | + | Ay | (ne + 1) | anes | 
+ | hy | Qna + 1) 1 engi | +... + 
+ 1 hy i(k — 1) ka + 1) 1 a ayn | - 


Let us assume that the estimate (4) is true for k = 1, 2, ...,(j — 1). Hence the above 


inequality yields, 
lh, | 
| A,||A, | (1+ oo) 


2n 


tL Aality | (1+ ThTY (it am) 
(j — In 





+ 


| hy | 
| jn(jnz + li amy | Stal +i al ae 


eet 
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which on simplification gives 
faa mL Eta ( L hy | 
7 [Ay (1+ : ) (1+ Fe) tet Data gs 1 P 
Lena jn (jna + 1) ° 


This shows that (4) is true for kK = j also. It is obvious from (5) that (4) is true for 
k = 1. Hence it is true for ¥ k > 1. 

Case (ii)—Leti A jn + 1 for j/>1. Then there exists a k such that (k — 1)n+1 
<i<kn+ 1. Here (5) reduces to 


i(ia+1—a)|a|Qlal +A lime + Ls ear 
+ | hy | (2nz+1) | Gonr | 
+ UAL IG 1) nee TI) Gyenae be 


Using (4) the above inequality becomes for (kK — 1)n +1<i<kn+ 1 








iia + 1—a)iari <i hy + yh) SEE 4 Lilly | Am 
n 
Lh |) | Ay | 
uae el Le Ost aren) | 
ill”, re aCe 


which on simplification gives 


hy 
og ee) 





Remark 2: When h is a general convex univalent function we can not ascertain 
anything about the sharpness of the coefficient inequality. However whenn = 2 
is 0 and h (E) is the half plane Re z > 0, we settlai<l, Yi>2, the sar 
estimate due to Sakaguchi*. When n = 2, « = | and h (E) is the half plane Re z > 0 


- 


we h i mat 
ave for ¥ k > 1, the sharp estimate | a, | < ka result of Das and Singh!. 


Theorem 4—We have C, (a;h) C C, (05h). 











PRooF : Let f(z) € C, (ash). Setti = ©) raphe 
we have a $n (Z) pep h ct ~ $n (2) 
$, (Z) 
a2 Uf" @))' + (1 — a) af" Gh SP CR Gem) 
oy eee SS 
azg,, (2) + (1 — a) dy (2) , (2) 
G)suanae 


(equation continued on p, 1121) 
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ES azp' (Zz) 
~?O+ G@+Fd-a* *O 


since f (z) € C,, (a;h). Here q (z) = h(z) (by our lemma). Again an application of 





lemma A yields p (z) = a on ~ h(z) which establishes the theorem. 


Now let us get an integral representation for functions belonging to the class 
Cn (%;A). 

Theorem 5—A function f (z) € C,, («;h) if and only if there exists a function 
G (z) in H(E) with G (0)=0 such that = 2 
with p (0) = | and p(z) 3 A(z) in Esuch that 





~ h(z)and an analytic function p (z) 


z 


f= sae | POEM 10m dr, if. x 0 


6 
and 
17) = peer Gh if a = 0, 
Proor : Let f(z) € C, (a;h). Then there exists a 


dE K, («;h) such that az (2f' (z))’ + CU — «) 2f’ (2), h (z) 
az oy, (z) + (1 — «) d, (z) 


where ¢, (z) is defined by (1). By our lemma ¢, € K, («;h); That is, 


az (2 ¢. (@)) +0 — a)z° (z) Pores 


azo (z) + (1 — «) 4, (z) 


Hence we can write 


(24, @) +0-826 ©  soiy 
1 5s, Se NG CE Pe Bee 





azg (z) + (1 — «) gn (Z) 


where 2G’ (2) -~ h (z) which on integration gives, 


G (z) 





az $' (z) + (1 —4) by (2) = G @) 
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and so, 





az (zf’ (2))'+ (1 — #) 2f'@) _ 
2Gf Cy p (2) = h(2) 


that is, ' | 
az (zf’ (z))’ + (1 — «)zf’ (z) = p (z) G(z). 


z(1/%)-2 





When « + 0 multiplying by and integrating we get 


z 


UG an [ro G (t) 101)? dt. 


0 


Conversely it is easily seen that if f(z) has the above integral representation then 
f(z) € C, (ah) 
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INVERSE AND REVERSIBILITY OF QUASIMATRIX TRANSFORMATION 
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(Received 6 August 1985) 


The object of the paper is to determine conditions under which the quasi 
method associated with an infinite matrix admits an inverse transformation. 


1. INTRODUCTION AND DEFINITION 


Vermes® pointed out that there is a close relationship between the summability 
properties of a matrix A = (a,x) regarded asa sequence-to-sequence transformation 
and those of its transpose A* = (a,,,) regarded as a series-to-series transformation. We 
call A* the quasi method associated with A. 


Given an infinite seriest 2 a, and the matrix A = (a,x), the quasi method A”* is 
to be ‘applicable’ if the A*-transformation 
Co 
b,, = 2% QAkn Ak ..(1) 
k=0 
exists foreach nm > 0. The series 2a, is said to be A* summable (or summable by 
the quasi-method associated with A) to the value s, if further 


xb, = S. td 4) 


When A is a normal matrix, then A* is in general arow infinite matrix and therefore 
it is relatively more difficult to work with. 


The generalised Norlund matrix (NV, p, 7) defined by 


Pack 4 (k a n) 
ank = * 

0 (kK >n) 
where r, = Po Jn +--+ + Pn Wo # O(n J YO), is normal provided that g, #0. In the 
case g, = 1 (all m) then the method reduces to the familiar Nérlund method (WN, p), 
see Hardy’ (p. 64). 

Kuttner® defined quasi-Cesaro as quasi-Hausdorff transformation. Thorpe’ 

has considered the transpose of the Nérlund matrix regarded as giving a series-to- 


series transformation in the case p, € DS, ie. 





+ Summation without limits means from 0 to o. 
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Pn > 0, a < ae <1 in 02 nee 
n n+] 


Mohapatro*’s has considered the transpose of the generalised Nérlund matrix (N, p. q) 
in the case p, © DZ, dn > 95 Ana > In- 


2. STATLMENT OF THEOREMS 


Generalising the results of Thorpe? and Mohapatro® we obtain the inverse of 
quasi-matrix transformation in Theorem | below whereas in Theorem 2 we determine 
conditions under which it is possible to recapture the original transformation from its 
inverse transformation. 


Before we state our theorems, we require the following notation : 


We write 4 € @ if A is normal and satisfies 


4 dn 
Opp 27. O(k-<ein);. Ae, SABE (0 kin I). 10943) 
Ansk+41 ank 


We write A € P* if AEC Pand 
Onk & Ankgi (0 Gk Cn — 1). ...(4) 


Theorem 1—Let A € @. Then the A* transformation (when applicable) has an 
inverse whose matrix is given by the transpose of the inverse of the A transformation; 
that is, if b, is given by (1) then 


Co 
a= > aj! be (8) 


k=n 


A matrix B, regarded as a series-to-series transformation, is called ‘reversible,’ if for 
any b = (b,) with = b, convergent, there is a ‘unique’ a = (a,) such that Ba = b. In 
the following theorem, we consider the reversibility of A* transformation, regarded as 
a series-to series transformation. 


Theorem 2—Let A € @*. Then A* transformation is reversible; more generally, 


if for any b, = o (1), then there exist a unique a = (a,) such that (1) holds. 
For the proof of the theorem, we require the following lemmas : 


Lemma | (a)—Let A € @. Then (see Peyerimhoff®, p. 33) 
CU} ox te, Ae oe 0; a, < 0(k > n); 


(b) Let4A € M*, Then (see Das?) 
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(ii) 33 | a,) | < © for each n > 0, 


k>n 


co 
(iii) SS a;} > 0 for each n > 0, 


k=n 


(iv) S la, |< 3 a,, (m > n). 


k=m+1 kan 


Lemma 2—Let A € #. Then for0 Gn < N < M, we have 








N 
= ay QGMo 
0< > a, ImMk S — < 


nn Ayo 
k=n 


PRooF ; Since by Lemma | (1), @,, > 0 and a,’ < 0 for k > 2, it follows that 


N M 
eT eR a > ~) ame = 0 
aan Mn = kn QAMk = Qin Mes 
k=n 


k=n 


by use of the identity (6) as M > n. But by (3) 


a fen - 
on eC. CL SRS rr). 


Ank-1 QAmk -j 


Hence it follows that 








A am | am oye amy 
a aMn Fame ” — . n . " Sad amo 
ut ann Ann aM-}; AMn-2 aMo 
cA 
< Bhi; Bee casei! Sal em ae 
Any any -1 nn =2 Ano 
_ @Mo 
no 


This completes the proof of the lemma. 


3. PRoor OF THEOREM | 
Substituting the value of b, as given in eqn. (1) and using the identity 


r 


La _ fle¢=n) 6 
S Den Chan Wa > n) Pal ) 


k=n 
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we have, for any integer V > n, 


N 


k=n ken rk 


N 2) 
=] 2 | 
Ay, OK = a, Ark Ar 


N r 
=) a | 
=> > ar > Qen ark + S ann > Ark Ar 


r=n k=n 


co N 


k=n 


; : ao 
=> an + ar Qin Ark. 


r=N+1 kon 


SCT} 


So it is obvious from (7) that (5) holds if and only if, for fixed 7, 


[ova] N 


-1 
Qn = > ar > a, ark > 0 as N > co, 


r2N+1 k=n 


co 
Since 5, = 2 Ako 4g converges by hypothesis, writing 


fore) 
Le 2 Ako Ak 
k=n 
we obtain 
WwW Ww 
a, — n n+l 
Ano 


Substituting this value of a, in (8), we have 


co N 


Were Wr 1 
Me aa a Pa a ark. 
. > Gro ae 


r=N+1 k™=n 
But for M > N 
M N 
Wr — W . 
SS “ Gr si » a, ark 
r=Nt+1 k=n 


...(8) 


.(9) 


...(10) 


(equation continued on p. 1127) 
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M N 
= a Qr- 
= } W, > an, ( me ut) 
Gro Gr .-150 
r=N+1 k=n : 
N N 
WN +1 ) -1 WM.) 4) 
—S— a. QnkE — a, AMk. ae 
any kn i amo kn Mk (11) 
k™=n k=n 


Since Ws is the remainder of a convergent series, it follows that w, = o (1) as k > c¢, 
Hence by Lemma 2, for each fixed n and N 


N 


Wu - ] 
aia Dy “in ame = 0(1)0(5° = 0 CD) 


k=n 


as M > oo, Also by identity (6), as N > n 


N 


a | 
S Gin aNk = 0. 


ken 
Hence it follows from (10) and (11) that 
=—>S Wr Dc Be (= _— Seat) ad 
(% Qr—150 
1=N+) k=n 


Hence in order that gv = o(1) as N — © for every (w,) such that w, = o (i) as 
r — oo, it is necessary and sufficient that 


2) S (2# = aa a;'| = 0 (1) ...(13) 


r@=N+1 k™n 





as N -— oo, for fixed n. 


By hypothesis (3) and Lemma | (a), the first term in the inner sigma in wy is 
negative and all other terms for k > are non-negative. Hence the modulus of the 


inner sigma in (13) is not greater than 
N 


a ar- = ark Gr —~15k —1 
v2 ( i aur) a= + > be ao ut) a (14) 
Gro ar _150 Arg Ar—150 


r=nt1 


and hence, it follows from (14) that 
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6a N 
ark ae ar 15k -—j | 
‘ > 5 (2 ) Ain a ) 


r=N+1 k=nt1 


By hypothesis (3), the first expression in (15) is 


= an : ay, 
a. ( Sa — lim ae) 
AN 50 moo 4mo 
a 
=] a. 
= mn o. = constant. 


Since the terms in the second expression in (15) are non-negative, we may invert the 
order of summation and the double sum becomes 





N oo 
> a SS ( - ) 
aro Ar -159 
k=ntl r=N+1 
N 
—] : ark a 
= S as ( line see 
aro ano 
k=ntl hats 
N 
] 
— aS ia ; Ss 
< He a, ANk (since 4,, <Ofork >n ) 
k=n+1 
=] 
ann G,,, 
ANo 


and this is bounded for fixed n. 


This proves (13) and completes the proof of the theorem. 


4. PROOF OF THEOREM 2 


Since 6, = 0 (1), it follows from Lemma | (b) (;; 
: th 
hence a, is defined by (5). (b) (ii) that (S) converges for each n; 


Now by identity (6) 
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N N 


co 
5 ay, Ay = } avy, ; a.) by 
vooae v=" k=y 
N N 
= Dom (St Dd ) ares 
vn key ke N +1 
co N 


oa 
= 
+ 
o 
>= 
g 
wat 
=) 
od 
< 


Hence 


N 
Hen = sae, a= pb. 
N->o von 


if and only, if, for fixed n, 


co N 
lim » by = ay, ay = 0. ...(16) 
esha k= N+1 v=n 


Since b, = o (1), to prove (16) it is enough to prove that, as N — oe, for fixed 7, 


an = >: > Av, aj = O(1). ...(17) 


But since by Lemma | (a), a, < 0(k > v), all the terms in (17) have the same sign 


and so we may omit the modulus sign. Thus 


fore) N 
-1 
aN — — avn ap 
k=@Nt+1 yon 
N co 
—1 
= QAvn —— any 
ven k=N+1 
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by (iii) and (iv) of Lemma 1 (b). 
This proves that (1) holds and the uniquenses of ‘a’ follows from Theorem e 
This completes the proof of Theorem 2. 


Remark : Theorem | is unchanged if we multiply each row of A by a non-zero 
constant and each column by a non-zero constant; in other words, we replace (4,«) in 
Theorem | by (An x @,x) where (A,) and (,,) are any sequences of non-zero constants. 


For suppose that A satisfies the conditions of Theorem 1. Applying that theorem 
with a, replaced by A, a, and 5, by b,/u, we see that if 


co 
b, = wn BU Ak Ak Ae 


k=n 


converges for all n, then 


co 


a = As > gh ay be 


k~n 


Thus Theorem | applies rather more generally, to any matrix of the form (A, ux @nx) 
where (4,).) satisfies the conditions in the form stated. 


ACKNOWLEDGEMENT 
We should like to express our thanks to Professor B. Kuttner for his valuable 


suggestions. 


REFERENCES 


_— 


. G. Das, Modified means (under communication). 

. G.H. Hardy, Divergent series, Oxford, 1949. 
B. Kuttner, J. Indian Math. Soc., 24 (1960), 319-41. 
P. Mahapatro, Pacific J. Math. 68 (1977), 177-95. = 
P. Mohaptro, J. Indian Math. Soc., 43 (1979), 259-69. 


A. Peyerimhoff, Lectures on summability, Lecture Notes in Mathematics No. 107, Springer- 
Verlag, Berlin. 


7. B. Thorpe, Ph. D. Thesis. Birmingham University, 1970. 
8. P. Vermes, Colloque sur la Théorie des Suites, Centre Belge de recherches Math., (1958), 60-86, 


NAnPv p 


Indian J. pure appl. Math., 16 (9) : 1131-1136, September 1986 


GENERALIZED STIELTJES TRANSFORM OF BANACH SPACE 
VALUED DISTRIBUTIONS 


A. K. Trwarli 
Department of Mathematics, Regional College of Engineering, Nagpur 440011 


(Received 1 August 1985) 


In this paper an extension of a generalized Stieltjes transform 


P= pet [ (4)® F(oe+15+0,- 5) 


f (t) dt; 
b=68+y+1,8 04> 0 
to Banach space valued distributions (generalized functions) is provided. An 


inversion formula for the above transform is also proved in distributional 
(Banach space valued) sense. 


1. INTRODUCTION 


The theory of distributions has been worked out in great depth by L. Schwartz. 
Sebastiaoe Silva' also has a theory for vector valued distributions. Zemanian?* has 
presented the theory of Banach space valued distributions. He has discussed the 
Laplace transform of Banach space valued distributions. Further he has used these 
concepts for applications in system theory and signals. 


The author*’® extended the generalized Stieltjes transform, namely 


re@+ynt+tiTe@+1)1 


ieee ep + 1) x 





co 
t 
| (<)PF (p+a+ HW Tass esi FE ea peg ae tan ed 2) 
0 
f (t) dt ati) 
to distributions. The complex inversion formula for the transform (1.1) namely 
c+wl r (d i 8 1) 
] s-68-— if 
= > Iii ae -____"_________ ¢-s M (s) ds 
ADH ges Be repays Ror er Lai (9) 


c—wt 
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where 


M (s) = hae F (x) dx 


b=PB4+nt+1d=b+a,s=ct+tiw,B >0,0>0 
was also extended to distributions by the author’. 


Motivated from the work of Zemanian*’* and Tiwari *’*’® we study generalized 
Stieltjes transform for Banach space valued distributions, in this paper. 


The notation and terminology of this paper will follow that of Zemanian”’*’? and 
Tiwari*’’® when U and V are two topological vector space, the symbol [U; ¥] denotes 
the linear space of all continuous linear mappings of U into V. A and e will usually 
denote Banach spaces. Throughout the sequel we use the following notation : 


R = set of all real numbers 
C = set of all complex numbers 
I = (0, oe). 


D (A) = space of smooth A-valued test functions with compact support defined 
on J. 


2. THE SPACES Sa (A) AND [Sa (A); B] 


Definition 2.1—The space Sa(A) is defined as the linear space of all smooth 
functions 9 (t) from J into A such that 


Past(?) = max sup || (1 + £)* t* oF (t)]l4 
O<k<! tEl 
<.00 5:1 <=» Opelge2in: 


The topology of Sa (A) is generated by sequence of seminorms {Pay}. - For 


A = Cwe write Sz (C) = Sy. It is easy to see that D(A) C Sa(A). If {a} be 
a monotonic sequence of real numbers such that a, — w+, where wis a real number or 
— oc, the space S (w; A) is defined as 


ioe) 
S(w;A)= U Sa (A). 
yl v 


For a Banach space B, any f € [Sq 


(A); B] will be called a Banach space valued 
distribution. When 4 = B = C, f becom 


es a scalar distribution. The simple or weaker 
topology for [Sz (A), B]is generated by the collection of seminorms {ye}, where o 
traverses S, (A) and 


ve (4) = I< f, e>I la. 
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We now come to the definition of the space [S. (A); B). 


If op € Sa, a € A, then pa = ag is the function from / into A that assigns to 
each t € / the valueg(t) a. Clearly p a © Sa (A). Wedenote by Sa @ A the sub- 


space of S. (A) consisting of elements of the form oa. Ifg € [Sa; C]anda € A, 
ga = ag is defined by 


<ga,9> = <g,9> 4,9 € Sa. 


Clearly ga € [Sa: A]. [Sa; CG] @ A denotes the subspace of [Sa; A] consisting of all 
linear combinations of elements of the form ga. Similarly when g € [Sa; [A; B]] and 
a € A, we define ga = ag € [S.; B]. Every g © [Sa (A); B] uniquely defines an 
f € [Sa; [A; B]] through the equation 

<f,o>a= <g,opa>, 9p E€ Sau aE A. 


See Zemanian’, p. 105. 


3. GENERALIZED STIELTJES TRANSFORM 

It is proved in Tiwari? that 
Tb T(B+1) 1 (4 i 
Td be 


F {5,8 +154; — ©) € Ss fora <1 


where 
b=B+yn+1,d=b)+ 4. 


Now let f € [Sa, A], we define the generalized Stieltjes transform F (x) of f by 


t 
FG) = =F): Ee) PGe) 1 (+ }" F(b,8 + iva; -4)>. 
rah) 


The left hand side F (x) is an A-valued function. 


Similarly if y € [Sa (A); B]. The generalized Stieltjes transform Y (x) of yis 
defined as the generalized Stieltjes transform of fy : 
)> 


...(3.2) 


x | 


jhe oy 1) 1 t \B a aoe 
Y (x) ay <fy (t), revert) (=) F (4,8 + Ay d; 


x x 


where 


fy € [Sa; [A; BI]. 
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: ae 
The above definition is meaningful because corresponding to y we have the unique fy 
as discussed in previous section. Note that Y (x) is an [A; B] valued function. 


4. AN INVERSION THEOREM 


We need the following there lemmas for the proof of inversion theorem. 


Lemma 4.1-—Let y € [Sa (A); B], fy be the corresponding member in [Sa; [A; B]] 


and 


6 (x, 7) = AM oe Gok F(,6-+15 4; — =). 


Then 
F x-* <y (x), 0 (x, > dx 
0 
od — 
=t <y(t), § 6(x,.c) <-* dx>. 
0 
Lemma 4.2—Let y € [Sa (A), B], 9 € D(A) and 
lo 2) 
P (s) = -J§ a(t) t-? de. 
0 
Then, for any two fixed real numbers r and c such that 0 < r < co,s =c + iw 


j <fy (t=), t91> P (s) dw 


= < y(s); f 2 Pls) dw. 


Lemma 4.3—If 9 € D(A) and a, c, r be real numbers with a<c. Then 


co ; ; « 
2”) ( - sin (; log dv 


u 
o ulog (=) 


converges in Sa (A) to 9 (u) as r > o0. 





Except for some obviuos changes proofs of the above lemmas are respectively 
similar to those of Lemmas 4. l, 4.2 and 4.3 in Tiwa 


ri*. We now come to inversion 
theorem. 


Theorem 4.1—If y € [D (A); B 
of Banach space valued distributio 
[D (A); B]. 


] and Y (x) be the generalized Stieltjes transform 
ny, thenin the sense of weak convergence in 
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c+ir 
; EE hehe AG es ona esas? ee 


r->o 


M(s) t-* ds, 9(t)> 


a 
my, o> 
where 


co 
M = { x*-! y (x) dx. 
0 


Proof : On substituting 


Tr(id+s—8-—1) LY, 
r(ib+s—-B—1)r(G@+1-—s)Ps 





and 
sS=c+ilw 
we have 


csir 


l T(d+s— B) 
Sowers ae DTG Ere oe ds, 9 (0)> 


c-ir 


iy 


eo = | OM (s) t-*, 9 (t)>. .. (4.1) 


—Ff 


The integral on w is an[A; 8] valued continuous function of t. By the continuity of 
the B-valued function Y (x) g(t) on the w — t plane and using change of order of 


integration (4.1) 


= a O M(s) [recoded 


lo 2) 
xs-1 y (x) dx | t-5 9(t) dt dw 
0 


} 
Pur 
eV" 
co) 
ares 


oO 
l eke Tib Ty (A) els foe" \0 
ae nee Mie) 
0 


(equation continued on p. 1136) 
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co 
aa ey ee t-5 0 (t) dt dw Bet 3" 
F (b,8 + 1s di zeae | 9 (t) 
] Ls oo 3 Pe (beat (=) 
- tl o<ve@, | x I'd x x 
Te 0 
co 
x F(b,8 + 1; d; — =) dx> | t~* p(t) dt dw <inte3) 
0 


by Lemma (4.1). 


Using the result of Tiwari‘ (p. 1052), viz. 


o[ x ean (z) F(%e+ 14-2) 
0 0 
x y(t) dt dx 


= fT ty (x) ds 
0 


we have 
r co 
(4.3) = ny | <y(t), v#-1> | t-5 » (t) dt dw 
=F 3 


r 


= <y (rt), a | 7s-1 | t~5 o(t) dt dw> (4.4) 


by Lemma (4.2). 


The proof now can be completed as in Tiwari* (p. 1053). 
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ON SINGLETONNESS OF UNIQUELY REMOTAL SETS 
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It is shown that a bounded subset G of a finite dimensional Hilbert space 


H, whose closure, G, admits a unique farthest point to its Chebyshev center 
is necessarily a singleton. The same result also holds when 4H is infinite 
dimensional and G is nearly compact. 


For a non-empty bounded subset G of a normed linear space X, the farthest dis- 
tance function from G is defined by 


t (x) = sup |ix — all, (x € X). 
a&G 


The subset G is said to be uniquely remotal (or to have the unique farthest point pro- 
perty) in X if for each x € X there exists a unique element g (x) € G such that 
t (x) = ||x — q(x)||. The following result is known for finite dimensional Hilbert 
spaces-. 

Proposition—Every uniquely remotal subset G of a finite dimensional Hilbert 
space is necessarily a singleton. 

We also know? that the unique farthest point property of a uniquely remotal 
subset G of a strictly convex normed. linear space is inherited by G, the norm closure of 
G. Knowing this, and the fact that Hilbert spaces are strictly convex, we can (equiva- 
lently) restate the above proposition as follows : 


Theorem —Let G be a subset of a finite dimensional Hilbert space H, such that 
G admits a unique farthest point to each x € H. Then G is necessarily a singleton. 


We now strengthen the above Theorem, as the following main result of this note. 
Recall that every bounded subset G of a Hilbert space has a unique Chebyshev center 


e (see Astaneh?, p. 1312'; i. ., there exists a unique c € H such that / (c) = bye 


(x). The non-negative real number t (c) is called the Chebyshev radius of G and is 
denoted by r (G). 


Department of Mathematics, University of British Columbia, Vancouver, 
Canada V6T 1Y 4. 


*Present Address : 
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Theorem 1—Let G be a bounded subset of a finite dimensional (real or complex) 
Hilbert space H. IfG admits a unique farthest point to the Chebyshev center c of G, 
then G is necessarily a singleton. 


Proor : Suppose G is nota singleton. Without any loss of generality we may 
assume that r(G) = | and c = 0 is the Chebyshev center of G. Let x, denote the uni- 
que farthest point for c = O in G. Then xp is the only point in G M U (A) with ||xoll 
= 1. (Here U (A) is the closed unit ballof H). Consider the sequence x, = p/n, 
(n > 2) in the open line segment ]0, x,[, and note that for each n 


1 <t(x,)? = sup ||x, — all. 
aGG 
(The strict inequality on the left is due to the uniqueness of the center 0 of G, which 
does not lie on JO, xo[). Therefore, for each n we can choose a, € G such that 
1 < |x, — yl? < ¢ (x,)’. ..(1) 
Clearly for each n we have |{a,|| < 1, and hence by (1), 


| ag lla,||? = Xn =a a,ll? ei IX nll? + 2.Re <Xny n> 
1 2 
= Ix, — al? — 4 +1(— ) Re <xy, a> 


l 2 
>1- nf + (—) Re < he; Gg 


Therefore, for each n the following inequality holds 


] 
2Re: = Xs, 0, = ea Ra 4 
Considering (1) and (2), for each n we have 


I|Xxo — aqll? = |] Xo — Xp Il? + x, — anil? + 2 Re <Xq — Xn, Xn -- Gn > 
alt 1 1 
—(1- 5) +i — al +2(1= =) 
‘ n 
RO ye 


1 2 7 
( n “iliac ian 2 (1 ut ~) n 


] 
~2(4 a 5) Re <x, a,> 


1 \2 
> (1 Se) ela I 
n n n (1- =) 


isn 
basis a n° 
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Therefore, for each n > 2 we get 


Ito—aull > a] 3) 


On the other hand since the farthest distance function x — ¢ (x) is continuous, and 
since x, 0 we have limit (x,) = ¢(0) = 1. From this and (1) it follows that 


x, — a@,|| > 1. Next from ||x, — a,|| > 1 and 
l 
Xp x a,|| oe Nn — Xn a a,| a lXnll < lla, || < l 
we deduce that lim |la,|| = 1. 


This means that {a,} is a maximizing sequence forO in G. We now Observe that, H 
being finite dimensional, {a,} has a cluster point 4 which of course lies inG. Since 


|@o|| = land |xXo — aol] = 5 (by (3)), it follows ‘that 0 has another farthest 
point a in G distinct from x,. This contradiction completes the proof of the theorem. 


We observe that the finite dimensionality assumption on H in the above proof is 
not used until the last paragraph, where we require a cluster point for the maximizing 
sequence {a,} for 0 in G. This observation enables us to deduce the following theorem 


in the infinite dimensional case, as well. 

Theorem 2—Let G be a bounded subset of an infinite dimensional (real or com- 
plex) Hilbert space H. Let G admit a unique farthest point to Chebyshev center c of 
G, and assume that every maximizing sequence for c in G has a cluster point (i. e. G is 
nearly compact). Then G is necessarily a singleton. 
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ON THE COEFFICIENTS OF POWERS OF A CLASS OF BAZILEVIC 
FUNCTIONS 
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Sharp estimates for coefficients of powers of a class of Bazilevic functions are 
given and the smallest subclass, yet known, of univalent functions, for which 
the powers of the Koebe function are not extremal, is introduced. 


INTRODUCTION 


Suppose that 0 < «a < oo. For functions f (z) =z + a,2° +a; 2° + 
analytic and univalent in the open unit disk U = {z: | z | <1} and for the Koebe ee 
ction k(z) = z (1 — z)-* we write 


po y - > a, («) 2" and / ‘oy = E b, (a) 2" (1) 


n=0 
n=0 


so that 


1 l -- = 
a, («) = Fae A i ABest | a, ASD a, Jana 


2a 
Meve 2(2+ a)(2+ io ss (2 + (m — 1) @) lad 
We consider the inequality 
| a, (@) | < B, (a). : et, 
If « = 1, this is the Bieberbach conjecture* which has recently been proved by 


De Branges* 24 for « = 2 it is the Littlewood-Paley' conjecture which was disproved 
by Fekete-Szegé5. 


In fact (3) is true for univalent functions if x < | and is false for « > 1 and 
n = 2, (See Hayman and Hummel‘), 


We may ask whether (3) is still true for subclasses of univalent functions. In this 
_Paper we introduce the smallest class, yet known, of univalent functions for which the 
Present Address : Department of Mathematics, University of Kentucky, Lexington, 
KY 40506-0027 U.S. A. 
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inequality (3) is false. Much of the work of this paper is drawn from the author’s 
D. Phil. thesis’. 


ror 8 > 0 let B(B), called Bazilevic of type 8, denote the class of functions f (z) 
analytic in U satisfying the differential equation 


2f’ (z) {Ff (z)}® 7? = ({g (z)}* p (z) ...(4) 
where g is a normalized starlike function so that 


g. (z) 7 
eae = q (2). wats) 


Here p (z) and q (z) are functions of positive real part with p (0) = q (0) = 1. It has 
been proved by many authors including Bazilevic? that the functions in B (8) are uni- 
valent. (See also Sheil-Small!'’?’). 


Note that B (1) is the class of normalised close-to-convex functions defined by Kaplan’. 


oo 
Theorem 1—Suppose that f(z) = z+ = a, 2" belongs to B(1). Then we have 
n=2 


the sharp bounds 





f 
ots fO<a<3 
! 
1 lla-3. 
leoaia ail) a 55 Ree 


Corollary 1\—For 3 < « <. eo, Theorem | shows that the inequality (3) does not 
hold for powers of functions in B (1). 


Baernstein! proved that iff (z) = z+ ... Is univalent, then for /y (z) 


raed . 
= {f(7pe4ozt+ & beet alibi) Se Kn-1|2, where K is an absolute constant. 
n=1 


t n—1+2/m 
Keogh and Miller? obtained tke estimate | Dinnti | < | ] for h,, (z) =z 


+ E Bog Or hae B(3)and0 <8 <1. Little is known about the co- 
n=1 


efficients of h,, (z) in 2 (m) when ™ is an integer and m > 2. We prove that 


Co 
Theorem 2—Suppose that /im( z) = z+ a Dans 2™"*! belongs to B (mm) where m 


is a positive integer. Then we have the sharp bounds 


ar Sane me ap 
} m 
| Brmar | < 4 llm — 3 


A Aaa i > 4, 
q 9m(m—1) eee 
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Corollary 2—For m = 2, Theorem 2 shows that Fekete-Szegd’s’ disproof of 
the Littlewood-Paley’® conjecture does not apply to the square-root transform of 
close-to-convex functions. 


Theorem 3—If mis an integer, m > 3, consider the meromorphic univalent 
function g (@) = £f (2")}m™ = E 4 cma 2" + Coma BM $n, 


where f(z) € B(1), z = 1/E and |&| > 1. Then we have the sharp bound 


llm + 3 
| eam < ~9m(m + 1) * 


If m = 1, the classical estimate | c:| <1 is sharp. The case m = 2 is more dif- 
ficult and we omit it. 


The estimate given in Theorem 3 corresponds to the estimate | cs,-1 | < m1 {1 + 2 
emt) m-1)) obtained by Fekete-Szegd* for the whole class of m-fold symmetric mero- 
morphic univalent functions. 


To prove our theorems we shall need the following lemma. 


oo 
Lemma |1—Let p (z) = 1 + 2 Pn z" be such that Re p(z)>0 in | z | < 1. Then 


S2+21/p,1? if > — 1/2. 





Ip. + a pt 

Proor : Write p (z) = {1 + @ (z)}/{1 — @ (z)} where w (Z=w,zt+wez-+ ... 
and | » (z)| <|z< 1. Equate the coefficients of z and z? to obtain P, = 2w, and 
P.= 2( We + wy ) By applying the fact that | w,.| <1 — |, | *we obtain | p, 


— 1/2 pp | <2— 1/2 ij pl? Nowfora > — 1/2 add(« + 4) | p,|*to the 


both sides of the latter inequality to obtain the lemma. 


Proof of Theorem 1—We see that the coefficients of functions f(z) in B(1) can 
be expressed by means of (4) and (5) in terms of the coefficients of two positive real 


art functi =. > « = = 
p ions p (z) = | api P, 2" and q (z) = 1 +2 gn 2” and so can the coeffici- 
ents a, («) of {z~1 f (z)}/™ defined by (1). 


Therefore we obtain a; («) = (pi + q,)/(2«) and 





Bee bys) 1 a+3 + 

a. = —_— — = teh 

= Tf gut gat st i a Pt BS ent 
...(6) 


It is easy to see that | a, (+) | & 2/x since | p, | € 2 and 14 |< 2. 
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For a, («) we obtain 




















l a+3 , l 3(« — 1) 
ela(o[< gat igra|+ sia *S Pe 
a + 3 
cs | Pa 
and so by Lemma | it follows that 
[—« a+ 3 da +3. 
a|a,(«)|< a | pi | + 6a Pilot eg IP ial 


We require the maximum of R (| p, | ) when 0 < | p, | € 2 and « is fixed and posi- 
tive. We note that R (0) = (7a + 3)/(6x) and R (2) = (a + 2)/a = aby (a). Thus the 


maximum of a, («) is greater than 5. («) if and only if R ( x) attains its maximum at PS 


in the interval (0.2). Now from 


d l—a a+ 3 
——__— = = 0, 
“capa 7 Prepare 2 Ak an may 








We obtain 
| po | = {2(@ + 3)}{3 (@ — I} andso R (| po | )=(I1«—3)/{9(a— D}. 


We observe that R(| p$ |) > max {R(0), R(2)} if3 < «<coandR (2)> R(0) 


if 0 < « < 3 since | P; | > 2 when 0 < « < 3. This yields the bounds in Theorem 


1. To show that they are sharp we write 9, = q. = P; = 2 and py = 2for0 <a 
< 3, p, = {2 (x + 3/3 («—1)} if « >3. This corresponds to p (z) = (1 + &)/(1—«) 
where w = z(pi + 2z)/(2 + piz) and q (z) = (1 + z)((I — z) and hence f (z) is con- 
structed from (4) and (5). 


To prove Theorem 2 we shall need the following well-known lemma. (See also 
Sheil-Small"?, p. 142). 


Lemma 2—Suppose that n is a positive integer. Then f(z) € B(B) if and only 
if { f(z)" € B (nb). 


Proof of Theorem 2—For hy, =A EVIE B(m) we obtain by Lemma 2 
that f(z) € B(1). Since for f(z) € B (1) both functions {z-1 f (ay and {z-"f 
(z™)}1/™ have the same coefficients which are defined by (1) form =a, we apply 
Theorem 1 to f(z) noting that bym4, = a,(m). Therefore we deduce Theorem a3 


Proof of Theorem 3—Since f(z) € B(\), its coefficients can be expressed in 
terms of the coefficients of two positive real part functions p (z) = | + pyz + --- and 
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qg(z) = 1+ qiz + ... and so can the coefficients of g (€) = { f(z™)}-’/". In fact with 
the notation of (1) we have c,,,-; = d,(— m). Therefore by applying (6) we obtain 


1 m — 3 m-+1 =, 
— MeCogm1 = 4 da + 3 Pat 24m qs 8m Py 


m— 3 
+ Tam Pa M 


Thus if m > 3 we obtain 








1 m— 3 m— 3 
m Com-1 < 6 I Qo | Bz 24m | q; | : -+- 12m | Pi qi | 
] 3 (mip 1) 5 
AEH etic mee 


Since 1/2 > 3 (m + 1)/8m, by Lemma | we deduce that 











l m — 3 im — 3 
m Com-1| S a moet | Py | F 2 6m P; | te ie GMEet 
m= RC [Py |). 
We observe that R (0) = (7m — 3)|(6m) and R (2) = (m — 2)/m. Also from 
d m+ | m— 3 
eae =— -———<<$_ ——_—_———-———- = 0 
d |p, | RAP, 1) 4m | pi | + 6m 


we obtain | p} | = (2m — 6)/(3(m + 1)}. Nowit is easy to see that R (| p, | ) has 


a maximum at | Py | because for m > 3, R( |p! 1) = (1Ilm + 3)) {9(m + 1)} 
= max {R (0), R (2)} and so the desired result follows. 
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FREE CONVECTION FLOW ON A NONISOTHERMAL FLAT PLATE AS 
INFLUENCED BY SUCTION OR BLOWING UNDER NONUNIFORM 
GRAVITY 


M. A. MANsour AND F. M. Hapy 
Mathematics Department, Faculty of Science, Assiut University, Assiut, Egypt 


(Received 15 October 1985) 


The effect of blowing and suction on free convection boundary layers deve- 
loped on a nonisothermal flat plate under nonuniform gravity are considered. 
A numerical solution is obtained for various values of poroisty parameter F,, 
and temperature distributions. The skin-friction and heat transfer rate on 
the wall at selected values of the distance from the leading edge and for some 
values of F,, are presented. 


1. INTRODUCTION 


Previous work on the effects of suction and blowing on free convection boundary 
layers has been confined almost entirely to the case of a heated vertical plate. Eichhorn! 
considered the power law variation in plate temperature and transpiration velocity 
which yields a similarity solution of the boundary-layer equations. Sparrow and Cess? 
discussed the case of constant plate temperature and transpiration velocity. They? 
obtained series expansions for temperature and velocity in powers of x1/4 where x is 
the distance measured from the leading edge. Merkin® extended the problem by obtain- 
ing asymptotic expansions, i. e. as x + oo, for temperature and velocity as influenced 
by both the blowing and the suction. The series expansions for small j given in Sparrow 
and Cess? were then joined to the asymptotic solution? by a numerical solution of the 


boundary-layer equations. 


In a recent paper, Clarke* expanded the analysis presented in Eichhorn! by obtain- 
ing further approximations to the solution of the full Navier-Stokes equations for large, 
but finite, Grashof number. In ref. [4] the density variations were considered, where 
as in Eiechhorn', Sparrow and Cess’, Merkin® the density variations were assumed im- 


portant only in producing the buoyancy force. 
Kao et al.> presented a solution for free convection along a vertical plate with 


arbitrary wall temperature variations in the absence of mass transfer at the wall. 
Aroesty et al.6 and Merkin’ considered the effect of blowing and suction through a 


body having a general shape. 


The purpose of this paper is to consider the effects of blowing and suction on the 
steady nonisothermal laminar free convection flows along (i) an infinite cold porous 
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plate rotating at speed w rads-' in a radial plane with its leading edge being ata dis- 
tance x, from the axis of rotation, and (ii) a finite hot porous plate of length xo, rotating 
at speed w rads~ ina radial plane about the line x = 0, assuming the gravity varies 
linearly with the distance x. The temperature distributions over the plate are also 
assumed to be polynomial of first and second degree, respectively. Numerical solutions 
were adopted and results for the skin-friction and heat transfer rate are presented for 
different poroisty parameter F,, (—0.43, —0.2. 0.0, 0.2 and 0.43) with Prandtl numbers 
Pr 0.7 and 1.0. 


The analogous problem of a free convection boundary layer on a vertical imper- 
meable flat pate under nonuniform gravity has been treated by, Pop and Na‘, and 
similar cases to those given in Pop and Na® are considered to solve the present problem. 


2, PROBLEM FORMULATION 


In the present analysis, an infinite cold plate rotating at speed  rads~’ in a radial 
plane with its leading edge beginning at a distance x» from the axis of rotations, anda 
finite hot plate of length x9, rotating at speed » rads~’ in a radial plane about the line 

x = O subject toa nonuniform gravity field g (x) are considered. The main effort in 
the present study is confined to the cases where the plate is porous and the plates tem- 
perature distributions 7, (x) have specific forms. 


Let the coordinates be chosen such that the x-coordinate is equal to the distance 
measured from the leading edge of the plate and the y-coordinate is the distance 
measured along the direction normal to the plate. The flow is assumed to be steady 
and incompressible. The boundary layer approximations are adopted along with the 
requirement that 2 (7, —7.)< 1 (large Taylor number). Thus the equations 
governing the conservation of mass, momentum and energy can be written in a dimen- 
sionless from as: 


du ov 
er Lee 0 weed) 
Ou Ou Ou . 
FS oT on = g (x) Sy (x) 6 + oy ...(2) 
66 d lin Sw (x) 06 1 a6 
uf peices a ]+°3 na 34 ay? i(3) 


with the boundary conditions 


=0,u=0,v= +y,(x), 9=1 
Pe sip ] ) 
Here u and v are the velocity components along the x, y coordinates, S\, (x) is the wall 


temperature, 6 the temperature, Pr the Prandtl number and v, (x) with minus sign 
corresponds to suction and with plus sign — to injection’ 
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Using the coordinate transformation 


Sy (x) 


e=x,y7 = “(axy oat 5) 


and introducing new dependent variables F and G such that 


y 
F(E1) = CHESTS, (xT? G (&4) = @. ...(6) 


These transformations reduce the boundary-layer equations to: 


FY + (3 + P(@®)) FF’ — 2[1 + P()) F? + 2(&)G 


=f (F = ate peal) (1) 
+ G’ + [3+ P@®] FG' — 4P(&) F'G = 48 (F a Go ) 
...(8) 
with the boundary conditions 
F' (E, 0) = + Fy = const., F’ (E, 0) = 0, G(E, 0) = 1 ] 
F~ (&, 0) = 0, G (&, cc) = 0 ...(9) 


where the primes denote differentiation with respect to y, and the function P (&) of the 
variable surface temperature is defined by 


dSw 
FG) = <5 See ...(10) 


The physical quantities of primary interest are the local skin-friction coefficient 
Cy and the local Nusselt number Nu which can be written respectively, in the form 


Nu = fee 


u 
Ou 


From the definition of wall skin-friction ty = »( “ey ) | yao and the local heat 


Cy = (11) 


transfer coefficient hy = gw/(Tw — Too), where gw = — k (@T/@y)y-o, Cy and Nu are 
expressed as follows 





3/4 
cp =4( SF) F&O, 


...(12) 





mu = -( & )" GG 0. 
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3. NUMERICAL SOLUTION 


In order to solve the differential equations (7) and (8), along with boundary 
condition given by equation (9), we assume that 


F= Fo(n) + € Fi(y) + & Fy (m) + ... ...(13) 
G = Go(n) + £G, (4) + &@ GL (4) +... .. (14) 


Up on substituting expressions (13) - (14) into the differential equations (7)-(9) 
and equating the coefficients of like powers of &, we obtain a set of ordinary differential 
equations. These are not reproduced here in the interest of conserving space. 


The governing equations for the velocity and temperature functions have been 
solved on the computer using the fourth-order, Runge-Kutta numerical integration 
procedure in conjunction with shooting techniques. 


The gravity g (&) is assumed to follow the form 
gs =-1+80GE cl #5(15) 


where the positive sign is taken for the cold plate and negative sign for the hot plate, 
respectively. 


4. RESULTS AND DISCUSSION 
4.1 Isothermal Porous Plate 


In this case S, = constant. The series solutions by Lienhard et al.° for the para- 
meters Cy and Nu can be written in the form . 


ep Aes APSSe 1m ara (16) 
n=0 nN 


Nu = ( ay" | - S (— 1)" ee é | (17) 


n=0 


where / = 0 for positive sign in equation (15) as/ = 1 for minus sign in the same 


equation. For comparison purpose, the calculates skin-friction and Nusselt number 
for cold and hot plates are given in Tables I and II res 
meter Fy, 


pectively at some porosity para- 
Figure | presents the dimensionless velocity F’ and the dimensionless temperature 
G for the isothermal plate (Sy = const.) at different values of Fy. Itis quite clear 
from this figure that an increase in Porosity parameter Fy leads to a decrease in F’, 
However, near the plate (n < 0.4) the temperature shows a decrease as F, increases. 
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Hot rotating porous plate Cold rotating porous plate 
9/9, =1-§ 





5 


Fic. 1. Temperature and velocity functions for the free convection along 
rotating isothermal porous plates with Pr = 0.7 and & = 0.2. 


Hot rotating Cold rotating 
porous plote porous plote 
G/on=1-3 gig, =1*8 


Hot rotating Cold rotating 
porous plate porous plote 


2.0 9! Qo =1-§ G!9,=1¢8 


FY,20,43 10 








0.0 
001 tt a0 OF tn OS 08 10 1.008 0.6 0.4 0.2 5 se 0.6 0:8 1.0 
$ 
Fic. 2. Comparison of F’ (&, 0)/F’” (0,0) Fic. 3. Comparison of G’ (€,0)/G’ (0,0) 
for the free convention along rotating for the free convention along rota- 
isothermal porous plates, Pr = 0.7. ting isothermal porous plates, Pr=0.7. 


Figures 2 and 3 refer to the relative changes of the local skin-friction fF" (& 0) 
F” (0,0) and the local heat transfer G’ (&, 0)/G’ (0,0) for representative values of & and 
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different values of ¥,,. For the cold plate the local skin-friction and the local heat trans- 
fer increase with the increase of &, but for the hot plate the effect of & is just the oppo- 
site. Also, the increase in the porosity parameter F,, leads to an increase in the local 
skin friction. The local heat transfer are identical. As is expected, g (&) hasa_signi- 
ficant effect on Cy and Nu as depicted by Figs. 2 and 3. 


4.2. Nonisothermal Porous Plate 
In this case, the class of wall temperature distributions defined, for instance as 
Sw(—&) =a, +a, +a, + ... ...(18) 


We shall assign values to the constants a;. Thus the problem is defined for the following 
cases : 


I Sy = 1 — 4& 
Tey ah bie hg epee (19) 


Figures 4 and 5 present the results for the skin-friction and heat transfer for the 
cases I —III evaluated at various & locations for different values of porosity parameter. 
Also, the profiles in the particular case of isothermal plates (Sy = const.) are also 


Cold fotating 
poras plote 
etaaet-t 


Het rototing Cold rotating 
porous plote porous picle 


Q'a,et-$ B'ggztef 


——Sy 2Const 





.0 0.0 
WO 08 060.4 02 0.0 02 04 0.6.08 10 


5 





Fic. 4. Skin-fraction results for nonisother- Fic. 5. Heat transfer resultsifor 
mal rotating porous plates. nonisothermal rotating porous 
plates, Pr = 0.7. 
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TABLE | 


Summary of numerical results for an isothermal cold rotating porous plate when P, = 0.7 


eed 


, C4/4 (Gre/4)>!4 Nu/ (Grz/4)'/4 
Fy = — 0.43 0,0 0.43 —0.43 0.0 0.43 
0.0 0.6630 0.6789 0.5579 0.2898 0.4995 0.8812 
0.0375 0.6788 0.6958 0.5772 0.2929 0.5043 0.8875 
0.1034 0.7064 0.7250 0.6111 0.2980 0.5124 0.8976 
0.1488 0.7251 0.7449 0.6344 0.3015 0.5178 0.9036 
0.2426 0.7632 0.7853 0.6826 0.3084 0.5283 0.9140 
0.3539 0.8074 0.8319 0.7397 0.3160 0.5398 0.9224 
0.4825 0.8571 0.8841 0.8056 0.3242 0.5518 0.9269 
0.5713 0.8905 0.9191 0.8510 0.3294 0.5593 0.9268 
0.6003 0.9013 0.9304 0.8659 0.3310 0.5616 0.9261 
0.6565 0.9219 0.9519 0.8946 0.3310 0.5616 0.9261 
0.7710 0.9631 0.9948 0.9530 0.3398 0.5737 0.8971 
0.9451 1.0234 1.0572 1.0417 0.3475 0.5834 0.8971 
ES 
TABLE II 
Summary of numerical results for an isothermal hot rotating porous plate when P, = 0.7 
Cy/4 (Gra/4)*"4 Nu|(Gr,/4))'4 

; Fy = — 0.43 0.0 0.43 —0.43 0.0 0.43 
pees ae eee ee 
0.0 0.6630 0.6789 0.5579 0.2898 0.4995 0.8812 
0.0375 0.6471 0.6619 0.5385 0.2867 0.9946 0.8743 
0.1034 * 0.6187 0.6316 0.5046 0.2811 0.4857 0.8612 
0.1488 0.5989 0.6105 0.4811 0.2772 0.4794 0.8513 
0.2426 0.5575 0.5662 0.4327 0.2687 0.4657 0.8286 
0.3539 0.5073 0.5123 0.3752 0.2581 0.4485 0.7978 
0.4825 0.4479 0.4484 0.3085 0.2452 0.4119 0.7256 
0.5713 0.4060 0.4032 0.2625 0.2359 0.4119 0.7256 
0.6003 0.3922 0.3883 0.2474 0.2328 0.4068 0.7148 
0.6565 0.3652 0.3591 0.2183 0.2328 0.4068 0.7148 
0.7710 0.3092 0.2985 0.1587 0.2136 0.3748 0.6453 
0.9451 0.2219 0.2037 0.0680 0.1928 0.3396 0.5644 


51 en te ti SSE 
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Mot rototing Cold rotding 
porous plote porous plate 


Hot rotating Cold rotating 


g/a,at-$ OlooeteS porous plote porous .plote 
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Fic. 6. Skin-friction results for nonisothermal Fic. 7. Heat transfer results for nonisothermal 
rotating porous plates case I. rotating porous plates case I. 


plotted at constant & distances in Figs. 4 and 5. For all the values of porosity para- 
meter, the skin-friction is greater in case I than in cases II and III of eqns. (19). For 
the cold plate the skin-friction increases with the increase €, but for the hot plate 
the effect of increasing — is to decrease the skin-friction. However for case III the 
skin-friction decreases with increase of § for both the cold and hot plates. 
also, the increase in Fy leads to an increase in the skin-friction. 
in Fig. 5 indicate that the heat transfer for all values of Fy is greater in case II than in 


cases I and III. Also, an increase in F, leads to an increase in the heat transfer for 
all the studied wall temperature. 3 


For case III 
The results shown 


Figures 6 and 7 show the calculated values of C s and Nu for case I at successive 
distances & for various porosity parameter with Pr = 0.7 and 1.0. The effect of P, is 


to decrease Cy and to increase Nu as Pr increase. The effect of Pr becomes predomi- 
nant at high values of F,,. 
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The present paper deals with the solution of the steady motion of a second- 
order fluid through the annulus between two coaxial right circular cylinders 
with suction and injection by finite element method. 


Galerkin’s method is used to reduce the linear boundary value problem into 
the corresponding set of algebraic equations. Gauss-Seidel iterative scheme 
is used to obtain the solution upto the accuracy of the order 1.0 X1078. The 
flow of the elastico-viscous fluid is taken as a particular case. The effects of 
second-order forces in the flow and that of suction and injection are discussed 
in detail and shown graphically. 


1. INTRODUCTION 


The study of flows in the presence of porous walls has many applications in the 
field of engineering, technology and biophysics. The viscous flow past a circular cylin- 
der and like problems are studied by many researchers'-4. Datta® has solved the pro- 
blem of steady motion of an idealized visco-elastic liquid through an annulus between 
coaxial circular cylinders and between two parallel plane boundaries, and flow of a 
non-Newtonian fluid through an annulus with porous walls. Sharma® has extended the 
problems considered in “’® for a second-order fluid. 


In the present paper we have investigated the Steady motion of an incompressible 
second-order fluid when it flows through the annulus between two coaxial right circular 
cylinders with suction and injection® by using finite element method. It is assumed 
that the rate of fluid injection at one of the surface is the same as that of withdrawal 
at the other, and the inner surface is moving with a constant velocity parallel to itself. 
The effects of the second-order forces in the flow and that of suction and injection 


tepresented by dimensionless parameters 7), T, ana R respectively have been discussed 
and shown graphically. 


2. FORMULATION OF THE PROBLEM 


Coleman and Noll’ have suggested following constitutive equation for an incom- 
pressible second-order fluid: 
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Cia a Pry + 2 by di; + 2 He C1; 4 Bs d* dx; (2.1) 
where 


di; — 3 (Vi, ; + V5) 
and 


Cij = 4 (4,,; + Ajst + 2 Ys Vins j) Azi2) 


p denotes the hydrostatic pressure and +,,; is the stress-tensor; while y, and a; denote 
velocity and acceleration vectors. pi, ». and ps are the coefficients of Newtonian-vis- 
cosity, elastico-viscosity and cross-viscosity?respectively. Equation (2.1) the momentum 
equation 


avi 
E ( ot + ,,v/) = ay Pleo) 
and the equation of continuity for steady motion 


y =0 eta) 


24 


where P is the fluid-density and a suffix following a comma denotes covariant differenti- 
ation, form the set of governing equations. 

Markovitz® has shown that 5.46 per cent solution of polyisobutylene in cetane at 
30°C behaves as a second-order fluid. The values of the material constants }1, », and 
j, are calculated by him are 18.5, —0.2 and 1.0 respectively (all expressed in C. G. S. 
system of units). On this basis the ratio of ». and »s has been taken as —0.2 for the 
purpose of numerical calculations. 

We work with cylindrical coordinates (r, 0, z) and choose the axis of boundaries 
involved to be along the axis of z. The velocity field for the problem can be chosen 
in the following form : 

u=u(r), v= 0,w = w(r). vein) 
In order that the velocity field (2.5) is compatible with continuity criterion we 


must have : 


d 
sets = 0 
dr (ur) 


which gives on integration 


ur = c, a constant. ...(2.6) 
Using eqns. (2.1) and (2.2) and velocity field (2.5) we have the physical compo- 


nents of stress tensor as follows : 
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du ve du \? (+ y"] 
Cir ee PF ht eee ty 2p ” det, oe +2 (St dr ) r dr 
du_\* dw \* 7) 
+e 4(4t) +(4) J 
u udu u* gre 
Tea = — DP a 2p aS + 2 }t2 | # a i I be Aus r2 6 (2.8) 
2 aw ...(2.9) 
to = — ptm 5) ( 
dw d?w du du dw 
Tr: = Py i oe Mo E dr2 T dr dr al ay 2 Hs dr’ “dr- 
mm 781!) 
Tro = Tos = 0 a2) 
and equations of motion for steady case are : 
Pu oe 4 a. Terr 49% bag ae ay Tee) tl eka? 
and 
dw d 1 
Pu etm mals 9 rail a(2.13) 


Let us suppose that z-axis is along the common-axis of the cylinders having radii 
aand b (a> b). If the inner cylinder moves with a constant velocity along the direc- 
tion of the axis and the outer cylinder is fixed, the boundary conditions for the 
problem are: 


ae US w= 0 at r= a 
am kW eet ae ...(2.14) 
From (2.6) and (2.14) 
aU = ur = Ub = a constant 


or, 


aU bU U 


Yu _—_—_—- TS —— ir 


r r y 
h On ee ; 
where 7 {| = a dimensionless distance. 


Equation (2.13) can be written as : 


dw d 
paU adr ed Tre] 
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or, 


paU 


Trz = ——_—. 
r 


(w + D) Byrne 


where D is aconstant. Eliminating tz, eqns. (2.10) and (2.15) give : 


d*w ages dw 
Bee ho, STi ea tak ONS: D) =0 (2.16) 


where 7; = p.U/u,a and T, = »3U/p,a are dimensionless elastico-viscous and cross- 
viscous parameters; R = UapP/p, is dimensionless suction parameter. 
The boundary conditions on w can be written as : 
w=0 when?7= 1, r(247) 
w= W whenyn= 6 


where s (= b/a) is ratio of the two radii. 


3. SOLUTION OF THE PROBLEM 


Ww 
W 


After dropping the asteriks, eqns. (2.16) and (2.17) can be written as : 


T; cat = (x? seth ere 223) ow R(w + D) =9 cl a Ba 
and 
w= 0 when 7 = |, ee eA) 
w= l when 7 = 6. 
Let the solution of eqn. (3.1) be 
w= z di (4) Wi ast(3.3) 


where ¢; (7) is linear shape function and defined as 
f 


ae when %j < 9 S 7141 


h 


1 (1) = 1 a 1 when qa <1 S 9 
| 
| 0 otherwise, 


where h is the step-length along %. 
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Putting the values of first and second derivatives of w and using eqn. (3.3), we get 
the error term as 


M M | ; 
E(y)-= T= $i (i) +E — (x acre eliey (n)w,— R 
{1 to 


M 
(2 $1 (qn) wi + D). ra) 
fo 
Applying Galerkin finite element method we get 


fE(). $41 =0,K=2 (1) N (3.5) 


ne, — T, — 27, 
ee aera =F | we. +( 22 — 7 Rh) WK 





h 2K-1 6 h 
2 
Mgt dy 2d 
+E OY Fd siege) — | 
he al ( 20K 41 ) a fet aeaaas 
K=2(1)N .-.(3.6) 
and boundary conditions (3.2) give 
w(1) = 0 wnen } = |, 
w(s) = ] when 7 = a. ia) 
Equation (3.6) in matrix form can be written as 
f Den sos ] f We = f RDh — a, w, } 
ey 5 
Bears ee RDh | 
| io | l 
| | | | | 
! : | | | 
| equal idan! | 
| an-2 5 cw-s Pi a | =e | 
{ Qn b } Wn J { RDh—cy_,wn3, J 
.-.(3.8) 
where 
45 geen 
ae ee Ahan ee ep 
Bes h 2ni-} 6 
b= mt — -—Rh 
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r. Nix, — Ty — 272 Rh 
a=— > +(—~——)- zim 1(I)N—1 
or, 
WEE 5 ( RDA Jel We 


OS) texts Wed) ) K=2(1)N. ...3.9) 


Computations are made for Wor (K = 2(1) N) by using the convergence cri- 
teria to achieve the desired accuracy (1.0 x 1078) 


Particular Cases : 


(i) by putting 7, = 0, we get the results for the elastico-viscous fluid. 
(ii) by putting the boundary conditions as 


Ww 


0 when 7 = 1 
w = 0 when 7 —co. 


The results are matched with Sharma et al.’. 


4. CONDITIONS FOR CONVERGENCE OF THE ITERATIVE PROCESS 
Equation (3.8) can be written as 


Ws = F, (We, Ws) 
wv = Fw-i1 (ww-1, W4i). 


S41) 


Then if wo ye AS ee be the approximate values of Wi, We.--, WN41 improved 


values are found by the following steps : 


0) (0) 
Hi? =F, (xf? 0) 


1 0) (0) 
WP =F (ww? ) 


(0) (0) ...(4.2) 
wi Fw (9, Pires ). ( 


Subtracting this system from eqn. (4.1), we have 


(0) (0) 
Wy — wv = F, (wi, Ws) — F, ( 6 » Ws ) 
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(0) (0) 
W3 —_— wo? = F, (W., W,) a F, ( w$ > W, ) 


\| 


Q) 
Les eee Wyn 


Fy (Wands WE ot EN (we, Wee (4.3) 


Now applying to the right-hand side of the system (4.3) the mean value of theorem 
we have 


: (0) (0) ' (0) OF, : (0) oF, 
F, (W,, s) -F( wf » Ws ) = (: — Wy, ) Ow, “ts (Wg se- Ws ) Cri 


where 

oF, OF [ w® 8 (Wj ww ) } wo + 0 (s - wo? )] 

Mea = hatte ee ee ae oe eee 0<e<l 
and 

oF, OF, [ wie = 7] (0, —¥ w ) b) wr + 7] (, + ww” 

Se Le ee Dye ae a 


OwW3 


In a similar manner we get 


Fz (Wo, W,) — Fa ( wo ee ) = (w, — wo ) or ~ (w. maT ) 





Ow, 4 
oF, 
ow, 
Fs (Ws, Ws) — Fs (no wm ) = (, — wi) ) 2S 
3 
a) 
ow; 


Fy_, (Wn-1, WN 41) — Fy-_, (ee : wo ) = (rns — Wena ) ru-s 
=I 


7 (was = Wee ) 
en there 2 
OWN 4) 
Substituting these expressions for the right-hand members of eqn. (4.3), it becomes 


Wy — Ww os (w, SH0 ) OF, i$ (x; Bens oF, 


ow, CN Be 
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maa (wa) Be alo ot) 
2 


4 OW, 
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@) (0) \ OFv-1 yo) OFw- 
WnN—- W = (ww. —_ YW ) tig 1) aed i w wr 0 1 : 
N 1 N-1 } Owy_, N41 — Wray he Ses 


(4.4) 


Adding these equations and considering only the absolute values of the several quan- 
tities, we have 



























































la ana we T | ws — wy fee n | < bw = w? | “a 
1 
+e. (0) | oe (0) oF, 
y |. We | aw, Ty Meade dws 
(0) oF, | ; (v) OF n_, 
Ss ion te aw, | +... bens Sa W arcs ae 
; eh OF y_, 1+(4,5) 
+ | WN41 Wy+4 | OWN a1 
Finally, as in Scarborough’, for convergence we have 
OF, OF; 
ew, OW3 ss I 
oF, oF, 
a ow, It 
| OFn-1 OFn-1 | | ...(4.6) 
| Own.) OWN 41 | ; 
Using eqns. (3.8) and (4.1), eqn (4.6) becomes : 
[bi <|e, | 
[b|<|alt+tla|,k=20)N—-2 
|b | < | an-1 | . .-(4.7) 


This is a condition for convergence of the solution. 


5. NUMERICAL VALUES AND DISCUSSION 


The numerical computations have been made for 7, = — 0.01, — 0.05, — 0.10; 
pee O04. = 018)5 = 51:2,-0.10, 0,15 0.20 and D = — 0.5, — 1.0, — 1.5, 0.05, 
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0.10, 0.15. Figure 1 depicts the behoviour of non dimensional component of velocity 
when 7, and R is fixed and the constant D is varying. It is noted that the behaviour 
of wis unaffected with D. The variation of w with suction parameter R when Lf, 16 
taken as fixed, is illustrated through Figs. 2 and 3 in case of suction and injection on 
the wall of outer cylinder. In first case, it is observed that the velocity component 
decreases as the suction parameter increases on the wall of the outer cylinder. In the 
case of injection, it decreases near the inner cylinder and increases near the outer cy- 
linder with an increase in R. Figure 4 shows the behaviour of w for different values 
of cross-viscous parameter T, when suction parameter is constant. It is found that the 
velocity component decreases with an increase in T, for suction as well as injection. 


__—_— - 
_-- 


--.- 





10 O4 os 
Le OCS i oes 
Fic. 1. Variation of Axial Velocity with D Fic. 2. Variation of Axial Velocity with 
for fixed 7; (= — 0.05) and R (= — 1.2). R (For suction case) for fixed 


T, (= —0.01) and D (= 0.05). 





Fic. 3. Variation of Axial Velocity with R Fic. 4. Variation of Axial Velocity with 
he injection case) for fixed T,; (For suction & injection cases 
T\(= — 0.01 )and D (= — 0.5). for fixed R & D 
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IMPLICIT NUMERICAL SOLUTION OF UNSTEADY 
EULER EQUATIONS FOR TRANSONIC FLOWS 


N. S. MADHAVAN AND VY. SWAMINATHAN 


(Received 11 November 1985) 


In the present paper, an implicit method described which makes use of the 
differential form of the unsteady Euler equations transformed to a rectan- 
gular computational domain. This technique entails a substantial reduction 
in computation time while, at the same time, retaining second-order accuracy. 
The method is applied to the case of transonic flow past a 6 per cent sym- 
metric circular arc airfoil and a computer program developed. The results 
from the software for surface pressure distribution on the airfoil have 
exhibited excellent agreement with the reported experimental results due to 
Knechtel for a 0.908 Mach flow, as also computed pressure values through 
an explicit scheme. 


NOMENCLATURE 
B = Jacobian of F, G 
c = Speed of sound 
D = Diagonal matrix 
e = Total energy per unit volume 
G = Vector flux (see eqn. (1)) 
I = Identity matrix , 
J = Jacobian of transformation 
n = Superscript representing advance of time 
p = Pressure 
= Used to define equation of body surface 
R’ = _ Seeegn. (4) 
= Diagonalisation matrix 
T = Temperature 


fs \== "Tinie 
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U = Dependent variable (p, pu, pr, e)? 
u, vy = Velocity components in Cartesian frame 
x,y = Cartesian coordinates 
A = _ Eigenvalue matrix 
P = Density 
—,» = Transformed spatial coordinates 
At, Ax = Elemental increase in ¢, x, ... 
Subscript 
i,j = i,j mesh pointsin &, 4 direction 
+ = Forward difference 


— = Backward difference 


Superscript 


— = Predicted value 


‘ = Transformed coordinate. 


1. INTRODUCTION 


Unsteady Euler equations have the advantage of being applicable in cases of 
supersonic, subsonic and mixed flows. Even in cases where other simplified flow models, 
such as steady or potential equations, are valid and could be solved without much dif- 
ficulty, the simulation of the actual flow from the initial conditions is of relevance not 
only because it is of considerable theorectical interest, but from the point of view of 
the insight it provides into the flow development. The solutions of these equations can 
be obtained by using the explicit time-marching method?’*’. Madhavan and Swami- 
nathan?*?, A major disadvantage of such methods, which are quite accurate, is, 
however, the stringent restriction on the time-step size and the consequent long com- 


putation time necessary to obtain convergence. 


MacCormack’s‘ implicit predictor-correct numerical scheme, which is second- 
order accurate, is easily applicable and involves a substantial reduction of computa- 
tion time by relaxing the restriction on the increment in time, may be made use of to 
overcome the above difficulty. The method involves only the inversion of bidiagonal 
matrices instead of tridiagonal ones as is normally required in implicit schemes. 
This technique is extended in this paper to the case of unsteady Euler equations trans- 
formed into a rectangular computational domain from the physical domain of any two- 


dimensional geometry. 
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The validation of the algorithm and the associated computer program was carried 
out by comparing the corresponding results obtained through the software with 
Knechtel’s experimental data for surface pressure distribution for transonic flow past 
a 6 per cent symmetric circular arc airfoil'. The computation time reduces to a small 
fraction (approximately 30%) of that necessary in cases where the explicit scheme is 
applied to the integral form of the unsteady Euler equations. 


2. GOVERNING EQUATIONS AND TRANSFORMATION 


The two-dimensional unsteady Euler equations in conservation form are given by 
ae op — 4 ——— = 0 eal) 


where 


——— 
a) im) 
cs = 


, 
| 
| 
| 
) 


soef 2) 


) 
| 
| 

G= | 
| 


Py? +p 
| 
(e+p)v J 
The numerical integration of these equations in the physical space is difficult 
owing to the fact that the grid will not be rectangular and consequently a transformation 
is made use of. It y = R(x) represents the equation of the body surface (eg. airfoil), 
the transformation § = x and y = y — R(x) transforms the computational domain 
into a rectangular one, while eqn. (1) assumes the form 
aU’ OF’ 
<r. SBE 
Considering the general transformation § = x and 1 = » (x, »); we have 
U' = UJ, F' = FiJ 


c= (Bee Bela $63. 


L 
c 
| 
{ 


eG’ 
Peappee ...Q3) 





In the present case, 
Ul se TF Foes 
G' = G — R’ (x) F. ..(4) 
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For a 6 per cent symmetric circular arc airfoil of chord length unity, the equation 
of the body surface, referred to rectangular coordinates with the x-axis along the 
chord length and origin at one end of the chord, is given by 


x*+ y? — x + 8.3033 y = 0 
and consequently the transformed coordinate 7 = x (x, y) is obtained as 


n=y — R(x) 


y — { — 8.3033 + +/8.3033 — 4x (x — 1)}/2. 


3. Tue ImMpecicir NUMERICAL SCHEME 


Differentiating eqn. (3) partially w. r. t. t, we obtain 


é oU 0 oU 7) aU 
or. i )= — GE (4 21 \= a Tm ) seo(9) 
where A and B are, respectively, the Jacobians of F and G’ w.r.t. U. Implicitly appro- 
ximating eqn. (5), we have 





eunti eu 
ie) kor a eu \"h F aU \nh 
At --=(4 at ) ~ £(85,) 
or 
Fy a aust 6U" 
(1 + At Eisoe + At font ) rym cael Tear re10} 
Defining 
aunt 
nt) —_ 
8U =, At 
and 
au" 
AUS ar At 
e sutt} — AU" (7) 
(7+ ar FEA + te B. | ntl _. AU™ . re 


Now an implicit predictor-corrector scheme may be defined for the numerical integra- 
tion of the transformed eqns. (3) as follows : 


Predictor : 





n Di, an Ds 
AU, =-“( aK Bi Gah Ay 


Gy ) 
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D. 
Ag 





D, n+l " 
(1-4 pal.) (7 ar R= 1B]. jog = aur, 


nm+r n > n+1 
Ce at UL 8 Le 





Corrector : 
= Fo Ee oe: OR Di eee 
AUT: —_ — ar( ay ta? ‘ pb An (Cri ) 
Kr eens r+ a: P= 131.) 30 = aus 
I+ At se | |.) re .) 87 is 
The Foe OR se OTe uz” | (8) 


The matrices | A | and | B | have’positive eigenvalues and are related to Jacobians 
Aand B. Let Se, Sn and their inverses denote the matrices that diagonalise A and B. 
Then, with perfect gas relations, A and B may be arranged as 











A= Sy’ AaSa ...(9) 
and 
B= Sy" AB Sn 
where 
a ea uB YB B } 
| Ch CC Ce 
| ue + aB C— up — vB B 
Sz = | 
ga I 
| p 0 5 0 | 
luc +o«B —c—uB — vB 8 | 
be up vB 8 
c C2 Cc? — oC 
Nx V yet, ny a Nx ‘ 
Se a Pas Pa, Pa, 


| 
| 
| 
| B, (ap —c¢ S-) a( - up) (Me - vB ) 8.8 
( 


8, («p aft ct) —8, (+4 ) ~Bi(™ c+ vB 8,8 


ied 


OS ee ed 
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f u 0 0 | 
| 0 ic 0 0 | 
Aa= | | 
0 0 u 0 | 
L 0 0 u—c J 
{ y 0 0 0 } 
| 0 chy 0 0 : 
As = | 
| 0 0 ay a CK, ) 
| 
a 0 0 0 4, — Cho 5 
2 2 LR 
a= Me + yy a2 = (a3 aaa ) 
ead 3 (u° am v*), B — ih y B, — 1]//2Pc. 5410) 
The matrices | A | and | B| are now defined by 
[Al = Sy Da Se 
and 
| Bi = S,* Ds Sa (iD 
where 
( AA, 0 0 0 : 
Ht ipoid. guadlier S07, bat ae 
| 
ra | 0 DagitiAcdet MOAN | 
| | 
0 0 0 Nag | 
f AB, 0 0 0 } 
| 0 ABs 0 0 | 
| 0 0 ABs 0 | 
l ® ; \By } ...(12) 





Z 


hea DE 
NA = Ma, = max{ | — 0.0 


I 


1 
AA, max {|utel— 3 gy 0.0 f 
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17 AG 
dag = max {iu —el- 5-4) .0.0 
1 Ay 
Ae = Aa, = max | | as | - 5 a 0.0 
1 A 
ray max { [a + cay |— 3-47 0.0 | 
1 A 
A pg = max | | a, — em | — sige 00H. 


(73) 


For regions of the flow in which At satisfies the explicit stability criterion, all A4 and 
Ag vanish and the set of difference equations (8) just reduce to corresponding explicit 
method. 


4. COMRUTER PROGRAM AND RESULTS 


A computer program was developed in FORTRAN IV language for carrying out 
the above computations and was made operational on the CDC CYBERI170/730 
Computer of VSSC. The experimental results used for comparison were those due 
to Knechtel for Mach 0.908 transonic flow over 6 per cent symmetric circular arc 
airfoil’. The upstream, downstream and top mesh boundary conditions were held 
fixed at uniform freestream flow in the expectation that the flow near the airfoil will 
converge before the disturbances reach those points which are held at 10 chord lengths 


away. The surface pressure on the airfoil was taken to be equal to that at the mesh 
point immediately above it. 


A 48 x 32 mesh system was used in the program. Along the x-direction, a two- 
mesh system of a stretched mesh (of 10 points) fore and aft of the airfoil and a uni- 
form mesh of 28’ points on the airfoil was made use of. Along the )-direction, a 
stretched mesh was applied. 


The data processing rate was about 8 secs. for iteration, i.e. 2.5 x 10-3 secs. 
per grid point including prediction and correction, that is, on an average, about 60 per 
cent more than that for the explicit scheme. Thus, when a CFL number of 5 was used, 
the implicit scheme took around one-third of the time for the explicit method. CFL 
numbers up to about 30 could be used without loss of accuracy. A maximum number of 
64 iterations were required to obtain convergence when the difference in consecutive 
flow parameter values becomes less than about 0.1 per cent (Shang and Hankey’) 
and this took around 4.0 minutes on the computer. The corresponding time for the 


explicit scheme applied to an integral from of the unsteady Euler equations for the 
same problem was approximately 11.0 min. 


Figure 3 presents the computed surface pressure values on the airfoil for Mach 
0.908 flow together with the corresponding experimental values due to Knechtel and, 
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Fic. 1. Mesh about typical circular ARC airfoil. 
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Fic. 2. Mesh system and boundary conditions in transformed computational plane. 
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Fic. 3. Comparison of pressure distribution for circular ARC airfoil (Moo = 0.908) 


on is quite good. The dotted line indicates 


as can be readily observed, the comparis 
ding cxplicit scheme. Figure 4 and 


surface pressure values computed by the correspon 
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Fic. 4. Pressure distribution for circular ARC airfoil (Mo. = 0.8). 
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Fic. 5. Pressure Distribution for circular ARC airfoil (Me = 1.2). 


5, respectively, exhibit the computed surface pressure for Mach 0.8 and Mach 1,2 
flows. 


5. CONCLUSION 


An implicit numerical method for solving unsteady Euler equations transformed 
to a rectangular computational domain is Presented in this paper. The method was 
applied for studying transonic flow past a 6 per cent symmetric circular are airfoil. 
and a computer program developed for the purpose. The method entails a substantial 
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reduction in computation time as compared to the corresponding explicit scheme, while 
retaining second-order accuracy. A comparison of the results for surface pressure dis- 
tribution obtained through the software for a 0.908 Mach transonic flow past circular 
arc airfoil with reported experimental values was quite good. Work is presently un- 
der way on solving Navier-Stokes equations using an analogous technique for the 
transonic shock wave boundary layer interaction problem. 
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